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ABSTRACT 

We present non-LTE time-dependent radiative-transfer simulations of supernova (SN) 
Ilb/Ib/Ic spectra and light curves, based on ~ 10^^ erg piston-driven ejecta, with and without 
^^Ni, produced from single and binary Wolf-Rayet (W-R) stars evolved at solar and sub-solar 
metallicities. Our bolometric light curves show a 10-day long post-breakout plateau with a 
luminosity of 1 -5 x 1 0^ L0 , visually brighter by ^ 1 mag than the progenitor W-R star. In our 
^^Ni-rich models, with '-^3 Mq ejecta masses, this plateau precedes a 20 to 30-day long re- 
brightening phase initiated by the outward-diffusing heat wave powered by radioactive decay 
at depth. A larger ejecta mass or a deeper ^^Ni location increases the heat diffusion time and 
acts to both delay and broaden the light-curve peak. Discriminating between the two effects 
requires spectroscopic modelling. In low ejecta-mass models with moderate mixing, 7-ray 
leakage starts as early as '^50 d after explosion and causes the nebular luminosity to steeply 
decline by ^0.02mag/d. Such signatures, which are observed in standard SNe Ilb/Ib/Ic, are 
consistent with low-mass progenitors derived from a binary -star population. We propose that 
the majority of stars with an initial mass <;2OM0 yield SNe II-P if "effectively" single, SNe 
Ilb/Ib/Ic if part of a close binary system, and SN-less black holes if more massive. Our ejecta, 
with outer hydrogen mass fractions as low as >0.01 and a total hydrogen mass of >0.001 Mq, 
yield the characteristic SN lib spectral morphology at early times. However at later times, 
~15d after the explosion, only Ha may remain as a weak absorption feature. Our binary 
models, characterised by helium surface mass fractions of ^0.85, systematically show He I 
lines during the post-breakout plateau, irrespective of the ^^Ni abundance. Synthetic spectra 
show a strong sensitivity to metallicity, which offers the possibility to constrain it directly 
from SN spectroscopic modelling. 

Key words: radiation hydrodynamics - stars: atmospheres - stars: supernovae - stars: evolu- 
tion 



1 INTRODUCTION 

Unlike numerous Type II SNe with a pr ogenitor identification on 
pre-explosion images (see lLeonardll2010l for a recent review), SNe 
Ilb/Ib/Ic can currently be constrained only through the analysis of 
their light. SNe Ilb/Ib/Ic events, understood as successful explo- 
sions following the core collapse of a H-poor/H-deficient progeni- 
tor star, represent a number of challenges for modem astrophysics. 
The traditional perspective was that, showing no hydrogen-line fea- 
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tures in their spectra, their progenitor massive stars had lost their 
hydrogen en velope prior t o explosion, in the form of a radiatively- 
driven wind JCastor et al.|[l975i) . During the star's evolution, mass 
loss progressively peels off the massive-star envelope leaving only 
a residual hydrogen shell. Eventually this is also peeled away leav- 
ing a shell whos e composition is dom inated by helium and ni- 
trogen (WN star; ICrowther et"aLlll995h . As evolution and peeling 
continue, the surfac e then becomes hel i um and carbon domi nated 
(yielding a WC star: lDessart et al.l2000l : ICrowther et al.l2002), and 
then carbon and oxygen dom i nated (yielding a WO star; Contil 
1 19761; iMaeder&Mevnelll 19941 ; iKingsburgh etai]l 19951) . Owing to 
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the stiff dependence of mass loss on luminosity/mass, low-mass 
massive stars lose little mass and die as Type II SNe. Consequently 
SNe Ilb/Ib/Ic would have to come from higher-mass (M <^ 30 Mq) 
massive stars. 

From the above simplistic argument, one can picture a se- 
quence of increasing main-sequence mass as we go through pro- 
genitors of SNe II-P. II-L, lib, l b, and Ic jHeger et alj booj : 
ICrowtheill2007L lOeorgv et al.ll2009l) . Recent revisions downwards 
of massi ve star steady-state mass-loss rates may challenge this 
scenario teouret et al]|2005l) . Although still quite speculative at 
present, transient, and perhaps recurrent, phases of intense mass 
loss (e.g., of the eruptive kind seen in 77 Car) may b e a viable al- 
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A second issue is the via- 
bility of the explosion mechanism for stars of increasing mass, 
since their higher-mass cores are increasingly more bound, expos- 
ing the newly-forme d SN shock to a tremendous accretion rate 
teurrows et al.ll2007h . Finally, we need to understand what differ- 
entiates the progenitors of SNe lib, lb and Ic, as well as the pro- 
genitors of Type Ic hypemovae which might be associated with 
a 7-ray-burst signal (fo r a discussion, see e.g. JWooslev & BloomI 
l2006l : lFrveretaLll2007h . 

Using radiati on-hydrodynamics simulations, 

lEnsman & WooslevI l ll988h were the first to perform a comprehen- 
sive study of hydrogen-less massive-star cores, representative of 
hydrogen-deficient W-R stars. Based on the available observations 
of SNe lb at that time, they concluded that ejecta from such W-R 
progenitor models could reproduce the general SN-Ib light curve 
morphology if powered by decay energy from ~O.1M0 of "'^Ni. 
However, the generally narrow peak of their light curves supported 
only 4-7 M0 ejecta, corresponding to 15-25 Mq progenitor 
main-sequence masses. Additional SN-Ib light curves acquired 
since the study of Ensman & Wo oslevI l ll988l) show a similar 
behaviour, and thus strengthen the notion that th e bulk of SNe lb 
are associa t ed with relatively low-mass ejecta jRichardson et al.l 
l2002ll2006l : lDrout et al.ll20Iol) . 

Since the scenario invoking higher-mass massive stars ap- 
pears not to be viable, binary-star evolution channe l is favored for 
the production of SN e Ib/c, as well as SNe lib jUtrobinI 
Wooslevet ai]|l994al : lYoung et alj 1 19951 : iBlinnikov etHl 
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A similar conclusion emerges fro m indepen- 
dent c onsiderations based on the observed SN Ib/c rate dSmith et al.l 
|20l3). 

In addition to single-star evoluti on models of W-R stars 
jLanger et al.ll994IWooslev et al.l2062l) . we now have physically- 
consi stent predictions for the binary-star counterpart (|Yoon et al.l 
20IO|: see also prelim inarv explorations of lWooslev et al.lll995l and 



Nomoto et al. I fl99l) . The critical ingredient of these scenarios is 



that a star may lose mass through mass transfer to a companion. 
Consequently, the mass-loss/luminosity scaling that inhibits mass 
loss in low-mass massive stars does not apply, allowing much lower 
mass progenitors to explode as "hydrogen-less" cores. 

The best observations and tailored analyses exist for only a 
few objects, usually the brightest or the weirdest, and thus not 
representative of the SN Ilb/Ib/Ic population. To cite a few, this 
includes the Type lib SN 1993J, associated with a small ejecta 
mass resulting from the expl osion of a rather low-mass progen- 
itor s t ar in a binary system dWooslev et al.l [l994bl : I Young et al.l 
I1995I: IBlinnikov et al.Hl998l). with a re sidual hydrogen envelope 
dSwartz et al.lll993al : iBaron et al.lll995l) . The Type Ic SN 19941 
was also the focus of numerous studies, suggesting a low-mass 



ejecta, potentially non-deficient in helium despite its classifi cation 
dWheeler et ai]ll994 iMiilard et al.l [19991 : ISauer et al.ll2006h . One 
exception is SN 2008D, which appears as a very standard Type lb 
SN but with a shock breakout detection, hence a very well-defined 
time of explosion. In this object, a short post-breakout plateau is 
observed before the SN re-brightens mer ely 4d after explosion 
and peaks 1 mag brighter about 20 d later ( Soderberg et al.ll2008l : 



IChevalier & Franssonil20"oi:lModiaz et al.l2009l) 

The helium abundance is thought to distinguish SNe lb and 
Ic, although the Hel-line signatures needed to make this as- 
sessment can be influenced by numerous complications. Due to 
the high excitation energy needed. He I-lines may not be ex- 
cited, and thus helium could be present but without associated 
spectral signatures. In addition, the presence or absence of He I 
lines is conditioned not just by the composition but also by non- 
thermal excitation/ionization processes bom out of 7- ray-emission 
from radioactive decay of ^^Ni and ^^Co isotopes (Luc 
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Kozma & Franssonll 1998 j lhT He I lines may therefore be affected 



by the efficiency of outward mixing of ^^Ni. Unfortunately, the 
mixing characteristics are adjusted for convenience rather than 
computed from first principles. Further, light-curve and spectral 
calculations are performed with separate codes that make distinct 
approximations and thus global physical consistency is lacking in 
cumnt modelling of SNe Ilb/Ib/Ic. 

Many SNe lb may have traces of hydrogen at the ir surface 
( iDeng et al.ll2000l : iBranch et al.ll2002l:|Parrent et al.ll2007h. and thi s 



may even apply to SNe Ic jjefferv et al.lll99ll:lBranch et al.ll2006l) 
While an interesting possibility, such inferences are always jeopar- 
dised by the weakness of the H I features at the heart of the debate 
(which suffer from l ine overlap with C II and Sill lines; see, e.g., 
iKetchum et alj2008l) and the shortcomings of the radiative-transfer 
approach used (e.g., UTE versus non-UTE, time-dependent effects 
on the ionisati on, ejecta structure and composition). As emphasised 
here (see also Ijames & Baronll2010h . obtaining early-time spectra 
and light curves of SNe Ilb/Ib/Ic is the only way to settle this and 
the He abundance issue. 

Here, we present non-UTE time-dependent calculations of SN 
ejecta based on piston-driven explosions of WN and WCAVO 
progenitor starsjj These SNe Ilb/Ib/Ic calculations have a higher 
level of consistency than earlier calculations — they combine 
stellar-evolutionary models of the progenitors, hydrodynamics cal- 
culations of the explosion (albeit artificial), and non-LTE time- 
dependent radiative transfer of the full ejecta yielding simultane- 
ously multi-color light curves and spectra. We discuss the gas and 
radiative properties of such SN ejecta and what determines their 
classification as SN lib, lb, or Ic. Our simulations ignore non- 
thermal electrons and do not include any mixing induced by the 
explosionl3 and thus set a lower threshold on the expected strength 
of Hi and Hel lines. More importantly, the calculations show that 
under certain conditions, which we detail. Hi and Hel lines can be 



^ W-R, WN, WC, and WO are spectroscopic designations. Here, for con- 
venience, we use the same designations for stars with the typical surface 
composition of a WN star etc. As some of the progenitors have very low fi- 
nal masses and large radii, but relatively low mass-loss rates, emission lines 
in their spectra may be unusually weak. 

^ In practice, the sharp ^''Ni-distribution left behind by the explosion is 
smoothed in our models by numerical diffusion as we remap the ejecta at 
the start of each step in our time sequences. This introduces a very moderate 
mixing, typically over a velocity-width of ~500km s^^. 
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seen even in the absence of any unstable nuclei and high-energy 
electrons. 

In Section [21 we give a brie f overview of the models selected 
from the comprehensive study of lYoon et al.l(l2010l) . including how 
we proceeded from such progenitors to make SN ejecta. As one 
aspect of this work is to understand how the progenitor properties 
effect the H and He lines observed in SN spectra, we discuss the 
difference in progenitor properties. We then present our re suits, dis- 
cussing simultaneously the properties of the gas (ejecta) and of the 
emergent radiation (spectra and light curves). The advantage of our 
approach is that it simultaneously yields synthetic light curves and 
spectra, allowing a direct assessment of spectra for a given light- 
curve phase. Our synthetic bolometric light curves and how they 
are affected by variations in photospheric conditions are described 
in Section [3] We then present our synthetic spectra, addressing 
in turn the potential signatures of hydrogen (Section 14.11 . helium 
(Section [4.2t . CNO elements (Section|43)> intermediate-mass ele- 
ments (IMEs; Section |4^ . and iron-group elements (Section [4.5l l. 
A crucial question which we address is when will the resulting SN 
be classified as Type lb, or Ic, or lib. Section |5] is devoted to a 
discussion of our results, which we confront to observations. This 
study is not a review so we merely select a few well-observed cases 
to illustrate our findings. In Section |6] we present our conclusions 
and future goals. 



2 PRE-SN EVOLUTION, EXPLOSION, AND MODEL 
SETUP 

The radiative-transfer calculations are based on SN ejecta that are 
produced in two steps. First, we simulate the evolution of a massive 
star from the main sequence until an advanced stage (either the end 
of neon-core burning or the formation of a degenerate iron core). 
Second, we use radiation hydrodynamics to simulate the explosion 
by driving a piston at the base of the SN progenitor envelope. 

For the stellar-ev olution model inpu ts, we use a sample of sim- 
ulations presented by lYoon et al] ( |2010|) . who focused on binary- 
star evolution for the production of W-R stars at solar and l/5th- 
solar metallicity and who used various primary/secondary masses 
and a range of orbital periods (typically on the order of a few 
days). In their study, the initial (final) mass of the primary was 
in the range 12-60 M0 (1.4-7.3 Mq), the final mass range reflect- 
ing the complicated mass-loss and mass-transfer history. These 
calculations were performed with an optimized nuclear network 
that included all elements up to y4 — 30, and iron. In a forth- 
coming study, a larger network containing all elements up to iron 
will be included. From this sample, we select a few W-R star 
models with a final mass in the range 3.79 to 7.3 M0, which we 
name here Bmil8mf3p79zl, Bmil8mf4p41zl, Bmi25mf5p09zl, 
Bmil8mf6p49zl, and Bmil8mf7p3z0p2. 

The adopted model nomenclature specifies the evolution chan- 
nel (B/S for binary-/single-star evolution), the main-sequence mass 
(which corresponds to X, given in M0, in miX; when present, "p" 
introduces the decimals), the mass at the end of the simulation 
(which corresponds to X, given in Mq, in mfX; this value should 
be very close to that at the time of core collapse in all cases), the en- 
vironmental metallicity (which corresponds to X, given in units of 
the solar value, inzX). Thus Bmil8mf3p79zl is abinary model for 
an 18 M0 main-sequence mass, with a final mass of 3.79 M©, and 
evolved at solar metallicity. We also complement this subset with 
two single-star models. These include Smi60mf7p08zl, a 20 M© 
helium-star model (60 Mq star on the main sequence) evolved at 



solar metallicity, and Smi25mfl8p3zOp05, a 25 M0 model which 
evolved chemica lly homogeneously at a metallicity of 0.05 Zq 
jYoon et al.ll2006l) . 

For historical reasons, this set of simulations is split into 
two distinct samples (see below). Models Bmil8mf3p79zl, 
Bmil8mf4p41zl, Bmi25 mf5p09zl were map ped as single stars 
and evolved with KEPLER jWeaver et al.lll978h until the formation 
of an iron core, before being exploded by means of a piston (using 
a '-^1.4Mq mass cut) to yield an asymptotic ejecta kinetic energy 
of 1.2 B (where 1 B = 1 Bethe = lO^'^ erg). The explosions produced 
0.171, 0.182, and 0.195 Mq of '^^Ni respectively, while the velocity 
of the ejecta shell that bounds 99% of the total ^'^Ni/^'^Co mass is 
2750, 2510, and 1250km s^^. A more detailed presentation cover- 
ing the post-helium burning evolution and explosion of these and 
other similar models is deferred to Woosley, Kasen, & Yoon (in 
preparation). 

The ID Lagrangian KEPLER simulations employ a grid of 
500-1000 points to cover the ejecta using a mass grid, thus resolv- 
ing and predicting sharp discontinuities in the composition distribu- 
tion. In contrast, our radiative-transfer simulations (performed with 
CMFGEN, see below), employ an optical-depth grid which is more 
relevant for radiation transport than a mass grid, but the drawback 
is that sharp composition boundaries cannot be resolved (at least 
without special procedures). As we step in time through our se- 
quence, the remapping procedure on such un-mixed models causes 
numerical diffusion of sharp composition discontinuities, which ef- 
fectively acts as mixing. We find that this smears the ®^Ni pro- 
file typically over 500km s~^. This also leads to a modest change 
in ^^Ni ejecta abundance compared to the KEPLER values, which 
are now 0.184, 0.170, and 0.237 Mq for models Bmil8mf3p79zl, 
Bmil8mf4p41zl, and Bmi25mf5p09zl (these cumulative isotopic 
yields are preserved at the 10% level throughout the time se- 
quence). 

The second set of models was used without furthe r evolution, 
remapped into VID l lLivnelll993l ; lDessart et al.ll2010bl lal). and ex- 
ploded, as in KEPLER, using a piston (with a 2.0 Mq mass cut) to 
yield an asymptotic ejecta kinetic energy of 1B0 Unstable-nuclei 
abundances were deliberately set to zero in this second set (we do 
not include any explosive nucleosynthesis for those VID simula- 
tions). Since these models were not evolved until iron-core col- 
lapse, the composition of the inner 2-3 Mq are not converged. Our 
discussion of such ^''Ni-deficient models will thus be limited to 
early times, when the photosphere lies well outside this inner re- 
gion potentially polluted by ^''Ni in realistic explosions. This dis- 
tinction between our two samples also serves to reveal the effect of 
decay heating on the light curves and spectra, which can only affect 
models Bmil8mf3p79zl, Bmil8mf4p41zl, and Bmi25mf5p09zl. 
It also permits a discussion of SNe lb light curves and spectra in 
which unstable nuclei are absent, which could occur, for example, 
with complete fallback of the inner ejecta. 

The splitting of the models into two groups primarily arose 
because the code of Yoon/Langer does not follow stellar evolution 
until the formation of a degenerate and collapsing iron core, and 

This slight difference of 0.2 B in ejecta-kinetic energy between these two 
sets of models has only a modest influence on our results. The difference 
in piston-mass-cut is of no relevance since we primarily focus on the early- 
time evolution which is sensitive only to the outer ejecta. The global prop- 
erties of our progenitor models (e.g., M,. R», etc.) would not be changed 
significantly if they had all been evolved until the fonnation of a degener- 
ate iron core, since the relevant nuclear timescales are shorter than, e.g., the 
Kelvin-Heknholtz timescale. 
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Table 1. Progenitor-model properties, including global characteristics and surface mass fractions of important species. We include the progenitor-envelope 
binding energy outside of 1 .4 (top three models) and 2.0 Mq (the rest of the models). Numbers in parenthesis refer to powers of ten. Ages generally refer to the 
evolutionary time since the main sequence. For model Smi60mf7p08zl, it is instead the time since the onset of core-helium burning. [See text for discussion.] 



Model 


Mi 


Z 


Ml 


R* 


Age 














^Mg ^Si 


^Fe 




[M0] [Zq] [Mq] 


[Rq] Myr 


[B] 
















Bmil8mf3p79zl 


18.0 


1.0 


3.79 


10.0 


9.8 


0.20 4.27(-2) 


9.38(-l) 2.02(-4) 


1.33(-2) 


3.76(-4) 


1.85(-3) 


7.25(-4) 7.34(-4) 


1.36(-3) 


Bmil8mf4p41zl 


18.0 


1.0 


4.41 


12.3 


9.5 


0.32 


1.25(-1) 


8.55(-l) 1.48(-4) 


1.33(-2) 4.53(-4) 


1.85(-3) 


7.25(-4) 7.34(-4) 


1.36(-3) 


Bmi25mf5p09zl 


25.0 


1.0 


5.09 


4.35 


7.5 


0.61 


6.33(-5) 


9.81(-1) 2.92(-4) 


1.33(-2) 


3.21(-4) 


1.85(-3) 


7.23(-4) 7.34(-4) 


1.36(-3) 


Bmi25mf6p49zl 


25.0 


1.0 


6.49 


3.64 


6.8 


0.58 


0.0 


9.81(-1) 4.23(-4) 


1.31 (-2) 


3. 11 (-4) 


1.69(-3) 


6.68(-4) 7.37(-4) 


1.36(-3) 


Smi60mt7p08zl 


60.0 


1.0 


7.08 


0.45 


0.5 


0.86 


0.0 


1.30(-1) 4.77(-l) 


0.0 


3.45(-l) 


3.42(-2) 


1.07 (-3) 8.25(-4) 


1.36(-3) 


Bmi25mf7p3zOp2 


25.0 


0.2 


7.30 


2.10 


7.3 


0.46 


0.0 


9.99(-l) 6.01(-5) 


2.6(-3) 


5.66(-5) 


3.23(-4) 


1.23(-4) 1.49(-4) 


2.72(-4) 


Smi25mfl8p3zOp05 


25.0 


0.05 


18.3 


1.14 


10.0 


1.49 


0.0 


5.34(-l) 1.72(-1) 


2.96(-3) 


2.17(-1) 5.82(-2) 


1.57(-2) 4.61(-5) 


6.80(-5) 



Table 2. Same as Table □ but now showing the explosion/ejecta properties computed by KEPLER and VlD, including the total ejecta yields for important 
species. The last column gives the velocity of the ejecta shell that bounds 99% of the total ^^Ni/^^Co mass. Because of the shaip composition boundaries 
in such unmixed ID models, these values are altered by up to 10% when remapping onto the lower-resolution grid in CMFGEN. [See Table[T]and the text for 
additional details.] 



Model 


-/\j?^rcmnant 






Mh 




Mc 




Mo 


A/nc A^Mg Ms; Mpe 








[Mq] 


[Mq] 


[B] 


[Mq] 


[Mq] 


[Mq] 


[Mq] 


[Mq] 


[Mq] [Mq] [Mq] [Mq] 


[Mq] 


km s~^ 


Bmil8mf3p79zl 


1.40 


2.39 


1.2 


1.94(-3) 


1.50(0) 


1.32(-1) 


1.19(-2) 


3.22(-l) 


9.24(-2) 4.68(-2) 7.94(-2) 8.94(-3) 


1.84(-1) 


2750 


Bmil8mf4p41zl 


1.50 


2.91 


1.2 


6.70(-3) 


1.64(0) 


1.63(-1) 


8.83(-3) 


5.36(-l) 


1.47(-1) 6.99(-2) 1.05(-1) 7.91(-3) 


1.70(-1) 


2510 


Bmi25mf5p09zl 


1.48 


3.61 


1.2 


1.39(-6) 


1.58(0) 


4.31(-1) 


8.56(-3) 


1.15(0) 


7.60(-2) 7.30(-2) 7.01(-2) 8.85(-3) 


2.37(-l) 


1250 


Bmi25mf6p49zl 


2.54 


3.95 


1.0 


0.0 


1.61(0) 


3.54(-l) 


8.08(-3) 


1.53(0) 


7.64(-l) 1.84(-1) 3.36(-2) 5.92(-3) 


0.0 


0.0 


Smi60mf7p08zl 


2.88 


4.20 


1.0 


0.0 


2.25(-l) 


1.19(0) 


0.0 


2.65(0) 


7.79(-l) 2.10(-1) 2.24(-2) 7.21(-3) 


0.0 


0.0 


Bmi25mf7 p3z0p2 


2.33 


4.97 


1.0 


0.0 


2.07(0) 


3.98(-l) 


1.21(-3) 


1.78(0) 


9.25(-l) !.22(-l) 4.19(-3) 1.40(-3) 


0.0 


0.0 


Smi25mfl8p3zOp05 


2.51 


15.79 


1.0 


0.0 


1.32(0) 


2.03(0) 


1.45(-3) 


9.11(0) 


2.83(0) 4.31(-1) 5.87(-l) 1.08(-3) 


0.0 


0.0 



because at the time VlD could not compute explosive nucleosyn- 
thesis. There is also a chronological reason for this split. The work 
presented in this paper includes computations started in 2008 (i.e., 
models Bmil8mf3p79zl, Bmil8mf4p41zl, and Bmi25mf5p09zl). 
Later, we realised that for a discussion of the early-time evolution 
of SNe Ilb/Ib/Ic, the actual core properties are not directly rele- 
vant. We thus augmente d our sample with additional models from 
lYoon et al. Il l2006l.l2010l) , evolved only until neon-core burning, and 
exploded them with V 1 D at this stage of evolution. Our choice of 
augmenting the sample was made after we had computed the light 
curve and spectra of models Bmil8mf3p79zl, Bmil8mf4p41zl, 
and Bmi25mf5p09zl. By including these ^^Ni-deficient models, 
we cover a wider range of properties, including higher-mass binary 
models, single-star models, solar- as well as sub-solar metallicity 
models, and a chemically-homogeneous model. 

In Tables [T] and |2l we summarise the properties of the progen- 
itors including the surface (outer ejecta) composition for the domi- 
nant species and the total ejecta yields, in solar masses, for impor- 
tant species. For each model we also quote the envelope binding en- 
ergy which is found to be signifi cantly smaller in bin ary-star than in 
single-star progenitors (see also lOessart et al.l2010bh - this is a par- 
ticularly important feature for the explosion mechanism (see Sec- 
tion|6j. The remnant mass Afrcmnant varies from small (^1.4]V[q) 
for the lower-mass models (evolved until iron-core collapse and 
exploded with KEPLER) up to 2.88 IVIq for the rest of the sample. 
The larger remnant masses may result, in part, from the flatter den- 
sity distribution of these model envelopes, a possible artefact of the 
truncated evolution which prevented the formation of a denser core 
suiTounded by a steeper density decline. Until we know how these 
stars explode the remnant-mass value we obtain may be inaccurate. 



At the start of the CMFGEN simulations, at 1-2 d after explosion, all 
these SN ejecta are in homologous expansion. 

Since the companion star in our selected binary systems 
may influence the SN radiation associated with the successful 
explosion of the primary star, we also give the binary system and 
secondary star properties. Binary-star models Bmil8mf3p79zl, 
Bmil8mf4p41zl, Bmi25mf5p09zl, Bmil8mf6p49zl, and 
Bmil8mf7p3zOp2 are characterised at the end of neon-core 
burning (of the primary star) by a systemic orbital period (separa- 
tion) of 29.7d (120.8Rq), 31.5d (125.6Rq), 21.3d (IIO.ORq), 
27.4d (127.9 Rq), and 13.0d (81.1 Rq). In the same order, the 
corresponding surface radius (stellar mass) of the secondary star 
is 10.9Rq (23.0]VIq), 9.3Rq (22.5]VIq), 26.0Rq (34.3]yiQ), 
II.SRq (31.0Mq), and 11.3 Rq (iSAMQ). 

Binary-star progenitor models tend to have larger radii than 
single-star ones, especially if they have retained a residual hy- 
drogen envelope. This concerns models Bmil8mf3p79zl and 
Bmil8mf4p41zl, which both contain a few O.OOl Mq of hydro- 
gen, with a typical surface mass fraction of 0.01-0.1. IVIodel 
Bmi25mf5p09zlalso has some hydrogen at the surface, but with 
a very low mass fraction of 6.33x10"^ and a cumulative mass 
of only 1.39x 10~® IVIq. After the explosion, such low-mass non- 
hydrogen deficient layers travel at speeds ^ 15000 km s^^ in this 
model, lie well above the photosphere past I d, and thus leave no 
spectral feature (Section 14.1) . The bulk of the envelope mass is 
lost through mass transfer, but the presence of hydrogen at the pre- 
SN stage is subsequently conditioned primarily by the stellar-wind 
mass-loss rate. In practice, the pre-SN model surface hydrogen 
abundance is higher for lower stellar-wind mass loss rates, which 
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Figure 1. Synthetic bolometric light cui'ves computed with our non-LTE time-dependent CMFGEN simulations based on binary- and single-star W-R progenitor 
models (see Table JTTll . A colour coding is used to differentiate the models, and dots indicate the actual post-explosion times at which the computations are 
performed (in a time sequence). Simulations including unstable nuclei (here exclusively produced through the decay of ^''Ni and ^''Co isotopes) brighten 
about a weak after the onset of a post-breakout plateau, while all other simulations show a precipitous fading as the photosphere recedes to the deeper helium- 
deficient ejecta layers of increasing mean-atomic weight and original binding energy. For reference, we give the bolometric magnitude on the right-hand-side 
ordinate axis. 



can result fro m the adopted ma ss-loss recipe, or from sub-solar 
metallicity (see lYoon et al.|[201ol for details). 

All progenitor models, except Smi25mfl8p3zOp05, have a 
rather small final mass, producing ejecta masses as low as 2.39 Mq 
in model Bmil8mf3p79zl. The low metallicity single star model, 
Smi25mf 18p3zOp05, produces an ejecta mass of over 15 M0 ejecta 
— over 3 times that of any other model. All models contain ~ 1 Mq 
of helium, with a surface mass fraction that is close to unity in our 
binary-star models, but substantially lower (0.13 and 0.54) in the 
two single star models. The model ejecta have a similar chemical 
stratification (reflecting the analogous progression through H, He, 
C etc. core burning), but they show quantitative difference in their 
envelope/ejecta yields. Consequently, we identify distinct spectro- 
scopic signatures associated with such abundance variations, which 
complement those associated with H and He. 

We evolve our hydrodynamical inputs usi ng the non-LTE 
time- dependent radiative-transfer code CMFGEN jHillier & Milled 
1 19981) . The approach and setup are analogous to that prese nted 
for the Type Il-peculiar SN 19 87A jPessart & Hillieill2010l) and 
for Type Il-Plateau (II-P) SNe jPessart & Hillieill201ll) ,"and are 
thus not repeated here. We assume the SN is free from external 
disturbances, such as interaction with the pre-SN wind or a com- 
panion star, or irradiation from the newly-born neutron star. Im- 
portantly, in ^^Ni-rich models, we assume a local deposition of 
radioactive-decay energy, treated as a pure heating source. Our 
treatment of decay energy as a pure heating source leads to an un- 
derestimate of the excitation and ionisa tion of the g as at the pho- 
tosphere as soon 7-rays may reach it (ILucvII1991|: Iswartzlll99lt 



(Kozma & Fransson"l992'; ' Swartz et al]|l993bl : iKozma & Franssorj 
il998a,b; Swartz et al. 1995). A corollary is that as long as there is 
no post-breakout re-brightening of our SN models, we are confident 
that the photosphere cannot be influenced by 7-rays from radioac- 
tive decay, and thus our computed spectra are not compromised by 
this approximation. Hence, we can discuss the early-phase spectro- 
scopic evolution of the SN when our non-LTE time-dependent ap- 
proach is physically accurate, starting as early as < I d after shock 
breakout. 

Our Monte-Carlo transport code predicts that 7-ray leak- 
age causes ~0.01% (10%) of the decay energy to escape at 
~40d (~85d) after explosion in models Bmil8mf3p79zl and 
Bmil8mf4p41zl. In model Bmi25mf5p09zl, this leakage causes 
~0.01% (10%) of the decay energy to escape at 75d (180 d) after 
explosion (Section l5.3t . Further, since LTE holds below the pho- 
tosphere, the bolometric luminosity is insensitive to how the 7-ray 
energy is downgraded. Thus, our synthetic bolometric light curves 
are strictly accurate over at least two months after explosion. 



3 BOLOMETRIC LIGHT CURVES AND EVOLUTION OF 
PHOTOSPHERIC CONDITIONS 

Starting our non-LTE time-dependent simulations at a post- 
breakout time of ^ 1 d, our synthetic bolometric luminosities ini- 
tially decline from the peak value reached at shock breakout. The 
early post-breakout fading stems from energy losses associated 
with radiation and expansion (the latter being particularly signifi- 
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Figure 2. Same as Fig.[T] but now showing the evolution of the radius (Rphot^ top left), the velocity (Vphot^ top right), and the temperature (Tp^ot; bottom 
left) at the photosphere. For reference, we also include the evolution of the total ejecta electron-scattering optical depth (Tojocta ; bottom right; the dotted black 
line indicates the 1/time^ evolution expected under constant ionization conditions). 



cant for compact progenitor stars), causing a tremendous cooling of 
the ejecta. In our simulations this rapid bolometric fading is largely 
missed — we catch mostly the end of it at 1-2 d. The luminos- 
ity then levels off at a value of 1-5 x 10^ L© and initiates a short 
plateau (Fig.lTJ. This plateau extends in all models over a number 
of days before the SN either brightens or fades. 

Brightening is initiated in our ''^Ni-rich models 
Bmil8mf3p79zl, Bmil8mf4p41zl, and Bmi25mf5p09zl as 
soon as the photosphere "feels" the influence of the heat wave 
generated from the decay of radioactive nuclei (primarily ^®Ni to 
^''Co) at greater depths. The re-brightening is delayed in model 
Bmi25mf5p09zl because its ejecta mass is larger and the ^^Ni- 
rich layers are located deeper (<1250km s~^) - both increase 
the diffusion time for the energy released at depth. In contrast, 
sample objects without ^''Ni eventually, and iiTevocably, fade at 
~IOd after exploding. This time is about the same as the time of 
re-brightening in the ^®Ni-rich models, and in some sense, is quite 
degenerate given the range of ejecta masses in our sample. The 
models that re-brighten reach a peak luminosity at about 30-50 d 
after shock breakout, before fading to a luminosity corresponding 
to the energy decay rate of ^®Co to ^®Fe (recall that we currently 
assume full-trapping of 7-rays in these CMFGEN calculations). 

This general light-curve morpho logy, discussed in a similar 
context bv lEnsman & WooslevI i ll988l) . is well known and gener- 
ally reflects the evolution of the photospheric properties, which 



we show in Figs. |2}[3] The modest increase in radius combined 
with the large decrease in temperature at the photosphere cause 
the luminosity to drop after breakout. The plateau arises when the 
temperature decrease at the photosphere slows down, which is as- 
sociated with the recombination of ejecta layers to their neutral 
state (pri marily He, but also CNO ele ments in such W-R progeni- 
tor stars). lEnsman & WooslevI d 1988b obtain a similar early post- 
breakout plateau in all their simulations. For reference, we also 
show in the bottom-right panel of Fig. |2] the time evolution of the 
ejecta electron-scattering optical depth, which tend to drop faster 
than 1/time^ due to changes in ionization. This figure also serves 
to identify the time at which each of the ejecta becomes optically 
thin. 

Free of any decay heating at early times, the evolution of the 
outer ejecta is essentially adiabatic (modulo radiative cooling near 
the photosphere), and thus ejecta internal energy is lost (and ion- 
ization state reduced) to pdV work primarily. What primarily de- 
termines the plateau brightness is therefore the amount of energy 
initially deposited by the shock and the size of the progenitor enve- 
lope. The former is modulated by the differential rate of ejecta ex- 
pansion between models. Given the uniform ejecta kinetic energy 
in our sample (I or 1.2 B), and the rough inverse correlation of size 
with ejecta mass for our limited set of models, models with smaller 
ejecta mass tend to have a larger plateau brightness (e.g., model 
Bmil8mf3p79zl has the largest kinetic energy per unit ejecta mass 
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Figure 3. Evolution of the composition (given as mass fractions) at the photosphere for models Bmil8mf3p79zl (top left; the time of peak brightness is 
~30d; model Binil8inf4p41zl has a very similar chemical stratification and is not shown), Bmi25mf5p09zl (top light; the time of peak brightness is ~50 d; 
the top axis is not shown since AMpiiot is a non-monotonic function of time in this model), Bmi25mf6p49zl (middle left), Sini60int7p08zl (middle right), 
Bmi25mf7p3z0p2 (bottom left), and Smi25mfl8p3z0p05 (bottoin right). The same colour coding is used in all panels of this figure to differentiate the species 
represented. Dots refer to the actual post-explosion times of the radiative-transfer computations (drawn on the curve corresponding to the helium mass fraction). 
Note that the timespan for ^^Ni-deficient simulations is limited to ~30d, while ''''^Ni-rich models (top row) are evolved all the way to the nebular phase. 



and the largest post-breakout plateau brightness). The latter is con- 
nected to the progenitor-envelope binding energy, which we show 
for a few models in Fig.|4] Objects that have brighter post-breakout 
plateaus are those with the least bound outer envelopes. Because of 
stellar evolution, these also tend to be more helium rich (greater en- 
richments in C/O testify for a greater proximity to the more tightly 
bound and hotter regions of the progenitor core where nuclear reac- 
tions take place). Consequently, ejecta produced both from a larger 
mass and a more compact progenitor star are the faintest of all dur- 
ing the short early-time plateau. 

The early-time plateau is not always strictly flat but instead 
may slant up or down. This seems to correlate directly to the 



depth variation of the ejecta helium abundance (and the associ- 
ated recombination energy) and of the corresponding progenitor- 
envelope binding energy. If the helium mass fraction is constant 
with depth (or if the progenitor-envelope binding energy varies lit- 
tle with depth), the SN luminosity is constant or increases, while if 
the helium mass fraction decreases with depth (or if the progenitor- 
envelope binding energy increases steeply with depth), the lumi- 
nosity ebbs. The main kinks in helium mass fraction shown in Fig.|3] 
thus correspond to points of inflection in the light curves of Fig.[T] 
In all ^"^Ni-deficient models, the precipitous fading at -^lOd coin- 
cides with the sudden decrease of the helium mass fraction at the 
photosphere, and its recession to the oxygen-rich, originally more 
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tightly-bound, ejecta layers. We find that the plateau is shorter for 
outer ejecta having a larger mean-atomic weight, because of the 
lower effective opacity and the greater original binding-energy (in 
absolute terms) of the corresponding layers. 

At the time of re-brightening in *^Ni-rich models, the heat 
wave from decay energy at depth causes the photospheric tem- 
perature to reverse its decrease and to rise. This rise causes a 
modest change to the ionisation, but helium remains essentially 
neutral at the photosphere (at later times, non-thermal ionisation, 
which we neglect, could alter this ionisation state). However, the 
recombination wave is slowed down and the photosphere recedes 
much more slowly in velocity space in models Bmil8mf3p79zl 
and Bmil8mf4p40zl (Fig.[2]l. During the rise to peak brightness, 
the photospheric velocity decreases by no more than a few percent. 
Because velocity relates to mass in SN ejecta, this implies that the 
photosphere probes the same region of the ejecta during that time. 
In model Bmi25mf5p09zl, the influence of the delayed heat wave 
is much stronger and causes the photosphere to move out in veloc- 
ity from 4000 to 5400 km s^^ (equivalent to an outward shift in 
mass of 0.6 Mq). 

At the onset of re-brightening the photosphere has receded 
through the very rich He-rich shell. While as much as < 50% of 
the ejecta mass may have passed through the photosphere, the 
photosphere is still well outside of the ^''Ni-rich regions. The re- 
brightening, and rise to peak, thus stems from the heat wave trig- 
gered by radioactive decay heating at greater depth. Because the 
7-rays are subject to a large optical depth at such early times, 
this diffusing heat wave will always reach the photosphere before 
direct heating by 7-rays can occur. More generally, diffusion of 
heat from the site of radioactive decay seems to be the cause of 
post-breakout re-brightening in all SN Il-pec/IIb/Ib/Ic progenitors. 
As the re-brigthening phase continues, the photosphere recedes 
through the He-C-0 region which is still He dominant, although 
in Bmi25mf5p09zl He, C, and O have comparable mass-fractions. 
At the time of peak brightness the photosphere is at the base of this 
He-C-O shell, and is entering the oxygen-rich core. 

Interestingly, the light curve of the ^^Ni-rich model 
Bmi25mf5p09zl qualitatively follows the early evolution of ^^Ni- 
deficient models, with the onset of a fading at 5-7 d as the photo- 
sphere leaves the helium-rich regions of the ejecta. However, this 
fading only lasts about a week before the heat wave reaches the 
photosphere and triggers a re-brightening. For yet a larger mass 
ejecta (with a similar ^^Ni distribution and mass), or for a deeper 
location of ''^Ni-rich material, this temporary fading could last 
longer. 

In ^'^Ni-rich models, the time to peak, the width of the peak, 
and the peak luminosity are conditioned by the ejecta mass, as 
well as the ^®Ni mass and its distribution. A larger ^''Ni mass en- 
hances the peak luminosity. Modulations of the ^^Ni distribution 
in mass/velocity space, as well ejecta mass (for a fixed kinetic en- 
ergy) alter the width of, and the rise time to, the peak, which in- 
crease with ^^Ni mass (amount of heating) or ejecta mass (diffusion 
time). The 2.39 Mq ejecta endowed with 0.184Mo of ^"Ni cause 
model Bmil8mf3p79zl to peak earlier (~30 d) at a larger value and 
faint faster after peak. The 3.61 M© ejecta endowed with 0.237 Mq 
of ^^Ni cause model Bmi25mf5p09zl to peak later (~50d) at a 
smaller value and fade slower after peak. Model Bmil8mf4p41zl 
behaves analogously to Bmil8mf3p79zl. 

As discussed above, the influence of decay energy (via the en- 
ergy diffusion wave) only kicks in at about 10 days in our ^^Ni-rich 
models (the photosphere hits the ''^Ni-rich layers not earlier than 
30d in our models; Fig.[3}- From this property, and prior to ~10d. 
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Figure 4. Illustration of the pre-SN envelope helium mass fraction (solid 
line) and binding energy (broken line), as a function of the Lagrangian mass 
(in units of the total ejecta mass) for our models Bmil8mf3p79zl (black), 
Bmi25mf6p49zl (turquoise), Smi60mf7p08zl (olive), Bmi25mf7p3zOp2 
(orange), and Smi25mfl8p3zOp05 (magenta). Models with brighter post- 
breakout plateau luminosities ai'e those whose progenitors have the least- 
bound outer envelopes, which also tend to be more helium-rich and stem 
from binary-star evolution. 



we can safely discuss the spectroscopic signatures in these simula- 
tions without any consideration of decay energy, gamma-rays, and 
non-thermal electrons, since these do not influence, either directly 
or indirectly, the emergent radiation. This also validates the use of 
progenitor models that were not evolved all the way to the forma- 
tion and collapse of the iron core — at early times only the outer 
ejecta are probed. 

As the spectral evolution is gradual over the first 10 d after ex- 
plosion, we discuss it at two representative times — 1.5 and 7.0 d. 
A montage of the spectra for all simulations is shown in Figs. [6]- 
[7] (where we plot x Fa for better visibility). From bottom to 
top, we stack the spectra in order of increasing ejecta mass. We 
have compiled numerous simulation results in these figures. Be- 
sides the full non-LTE time-dependent spectrum, we also show as 
coloured lines the spectrum obtained when the bound-bound tran- 
sitions of a selected species are omitted (for example, the red line 
shows the spectrum resulting from the neglect of H I-bound-bound 
transitions). This does not give the line flux associated with a given 
species, but clearly marks the spectral regions that it affects. Fur- 
thermore, below each spectrum, we draw tick marks for all lines 
with an equivalent widths EW greater than 10 A, with a thickness 
that scales with log(EW). For this figure, we address the contribu- 
tions from the following important species: H I, He I, C I, N I, O I, 
Na I, Mg II, Ca II, Si II, Ti II, and Fe II. 

At ~1.5 d after explosion (Fig.|6}, our synthetic spectra show 
a similar slope, with the peak of the spectral energy distribu- 
tion (SED) at ~5000A. Line features are more numerous and 
pronounced in the blue part of the spectrum, although a few 
broad features are also present in the red part of the spectrum. 
As indicated in this montage, we identify lines of Hi (models 
Bmil8mf3p79zl and Bmil8mf4p41zl), of Hel (all models except 
Smi60mf7p08zl and Smi25mfl8p3z0p05), and of once-ionised 
CNO elements. The standard CallH&K and 8500 A-multiplet 
lines are present, although with a strength that varies depending 
on the Ca ionisation state. A forest of lines from Ti II and Fe II 
are present in all models apart from the low-metallicity models 
Bmi25mf7p3z0p2 and Smi25mfl8p3zOp05. Nearly all the lines 
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Figure 5. Same as Fig.|2] but now showing the evolution of the mass above 
the photosphere (AA/pij^t)- Notice the sizeable reversal of the photosphere 
trajectory in velocity/mass for model Bmi25mf5p09zl. 

we predict in model Bmi25mi7p3z0p2 are from He I. In contrast, 
model Bmi25mf5p09zl essentially shows just Fell lines. At 7.0 d 
after explosion (Fig. |7}, the SEDs are significantly reddened and 
line blanketing effects are stronger, primarily because of Fe II. 



4 SPECTRAL PROPERTIES DURING THE 
POST-BREAKOUT PLATEAU 

Overall, the different spectral properties are conditioned by the dif- 
ferences in composition between models, which reflect the differ- 
ing evolutionary paths and environmental metallicity. The compo- 
sition is also important because ejecta elements have a range of ion- 
isation potentials. It is about 1 1.3-13.6 eV for H I, C I, and O I, but 
it is 24.6 eV for He I, thus about twice as high. Due to the high ion- 
isation potential of He it may not contribute to the spectrum, and 
thus its presence can be hidden. Furthermore, the ionisation state 
of the photosphere is time dependent and conditioned by complex 
collisional/radiative processes, requiring detailed radiative-transfer 
simulations. Below, we present our results for lines associated with 
the main elements of interest: H, He, CNO, intermediate-mass ele- 
ments (IMEs), and iron-group elements. 

4.1 Hydrogen 

In our model set, only models Bmil8mf3p79zl, BmiI8mf4p4IzI, 
and Bmi25mf5p09zl contain hydrogen, with a maximum mass 
fraction of 0.043, 0.125, and 6.3x10^'', and cumulative masses 
of 1.94xl0"^ 6.7xl0"^ and 1.39xlO"^M0. This hydrogen is 
typically present above ~10000km s~^. In the corresponding re- 
gions, the composition is 85-98% helium, with an integrated mass 
of ~I.5Mq. In model Bmi25mf5p09zl, the hydrogen mass frac- 
tion of ^IO~*Mq is too small to produce any feature, but in the 
other two models, synthetic spectra clearly show Ha. Other Balmer 
lines are weak and overlap with metal lines which can compro- 
mise their definite identification (for example, H/3 merely broadens 
a broad feature due to iron-group elements). At 7.0 d after explo- 
sion, Hq becomes the second strongest/broadest line in the opti- 
cal spectra of models Bmil8mf3p79zl and Bmil8mf4p41zl. It is 
stronger and broader in the latter model due to its larger hydro- 
gen mass fraction. Interestingly, our early-time synthetic spectra 



for models BmiI8mf3p79zI/Bmil8mf4p41zl look very similar to 
those of Bmi25mf5p09zl if we were to simply remove the Ha line 
from the former models. This suggests that such low amounts of 
hydrogen do not influence sizeably the radiative-transfer solution, 
but can alter the SN classification. 

The presence of H I lines (and primarily Hq) calls for a Type 
II classification. The 7-day spectrum does not look quite like that 
of a SN Il-Plateau at the same epoch jPessart et al.ll2008h . and the 
corresponding model light curve even less so. The subsequent evo- 
lution exacerbates the differences as the strength of the H I lines 
rapidly decreases as the photosphere recedes to H-deficient layers 
of the ejecta (Fig. [3). Both H/? and Ha vanish at ~ 10 d after the ex- 
plosion. Prior to that, Hq shows a broad flat-top profile with a weak 
absorption feature at large velocity (Fig. |8). With this morphol- 
ogy, the line is often referred to as being "detached". This stems 
here from chemical stratification, with hydrogen being absent in all 
ejecta mass shells with a velocity ;^I0000km s~^. In the present 
context, and if helium lines strengthened at later times, this event 
would accurately be classified as a Type IlbQ 

The evolution of Balmer lines reflects in part the ionisation 
conditions at and above the photosphere (Fig. [3}. As time pro- 
gresses, the regions in the vicinity of the photosphere recom- 
bine and cool, photoionisation/recombination rates decrease, so 
that the Balmer lines eventually weaken. This is mitigated at 
large radii/velocities by an ionisation freeze-out which stems from 
time- dependent effects dUtrobin & Chugaill2005l;lDessart & Hillieil 
I2OO8I) . The arrival of the heat wave at the photosphere at ~8 d after 
explosion raises the hydrogen ionisation and causes the hydrogen- 
recombination front to travel outwards in velocity /mass. Because 
hydrogen is under-abundant in those ejecta regions, it does not 
reverse the recession of the photosphere to deeper ejecta layers, 
since helium, the dominant species there, continues to recombine. 
Spectroscopically, this ionisation shift causes the H/3 absorption 
to plateau, while in Ha, it causes a kink in the P-Cygni absorp- 
tion. The Ha-line survival time, strength, and width, could easily 
be modulated by varying, for example, the progenitor radius, forc- 
ing hydrogen to recombine earlier/later, or by mixing hydrogen to 
deeper layers (this coud result from stellar evolution or from fluid 
instabilities triggered after shock passage). 

From these simulations, which are based on realistic binary- 
star evolution models, we predict that H I lines are produced for 
very low hydrogen abundances, corresponding here to ^O.OI mass 
fraction (if <I0^'*, no Hi line is obtained) and >0.001 M0 cumu- 
lative mass (these values appear quite standard for the correspond- 
ing progenitor regions, which are at the base of the hydrogen en- 
velope and thus show significant mixing of hydrogen and helium). 
The Hi lines are predicted for up to ~10d in the absence of non- 
thermal processes and more generally without any need for ^^Ni. 
Hq, the strongest Balmer line we predict, is strong at early times 
but considerably weakens as time goes on, so that its unambiguous 
identification can only be done at early times. At later times, it may 
still be present but this will be conditioned by multiple, compli- 
cated, and sometimes largely unconstrained effects. These include 
mixing, ^''Ni production and distribution, ejecta mass, chemical- 
stratification etc... Further, its identification will be compromised 
by overlap with numerous and potentially stronger lines. Obvi- 



* SN classification is, however, subject to errors. For example, SN 2008ax 
was origin ally classified by Blond in et al. 1 2008) as a Type Il-pec and later 
revised by IChomock et al.l j20O8l) . based on additional observations, as a 
type lib. 
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Figure 6. Scaled synthetic spectra (black) for each of our simulations at ~ 1 .5 d after explosion (we show the quantity X _F>, to better reveal the weaker 
features at longer wavlengths). At each epoch, we oveiplot the continuum-only synthetic flux (black dotted line; the same vertical scaling is applied), as well 
as the synthetic spectra obtained when bound-bound transitions of a given ion are omitted (see colour-coding at top). This highlights the relative importance 
of lines, and their associated ions, in different spectral regions. Coloured tickmarks indicate the ion producing a line with a Sobolev equivalent width of 10 A 
(in absolute value), and with a thickness that reflects the magnitude of this equivalent width on a logarithmic scale. 
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© 2010 RAS, MNRAS 000.[Tll23l 



12 Luc Dessart et al. 




Time Since Brealtout [d] 
1.5 4.6 7.6 10.6 13.7 



-0.6 b 



-20000 -10000 10000 20000 

Doppler Velocity [km s~'] 




Time Since Brealiout [d] 
1.5 4.6 7.6 10.6 13.7 



-20000 -10000 10000 20000 

Doppler Velocity [km s~'] 

Time Since Breakout [d] 
15 22 




10000 



30000 
[km s"' 



40000 



50000 



20000 
Velocity 

Figure 8. Top and Middle: Evolution of H/3 (top) and Ho (middle) line 
fluxes versus Doppler velocity from 1.5 to 13. 7d after shock breakout for 
inodel Bmil8mf3p79zl. More specifically, the quantity plotted is obtained 
by subtracting from the total flux (F\) the flux obtained by ignoring all 
bound-bound transitions due to H I (FJ^ ) , and then normalizing by the con- 
tinuum flux (Fc, i.e., the model flux obtained by accounting only for con- 
tinuum processes). This approach illustrates the H I-line flux contiibution to 
the total flux and its evolution with time. Bottom: Illustration of the mean 
hydrogen ionisation level for model Bmil8mf3p79zl. The coloured dots 
refer to this ionisation level at the photosphere and its evolution with time. 
Notice the re-ionization of hydrogen above the photosphere after >7 d in 
model Bmil8mf3p79zl caused by outward-diffusing decay energy released 
at depth, and which causes a kink in the Hq absorption profile (middle 
panel). 



ously, obtaining very early-time spectra alleviates all these issues 
and can set a firm constraint on the properties of the progenitor 
star. 



4.2 Helium 

The production of He I lines from ejecta arising from W-R-star 
explosions has b een routinely associa t ed with non-thermal exci- 
tation processes jHarkness et al.l 19871: Lucvlll991 : Swartj T99ll; 
Kozma & Fransson'l992':' Swartz et al.lll993bl ; Kozma & Franssoiil 
T998a,b; Swartz et al. .199^ This has been motivated by the in- 
ability of producing them by other physical means. However, these 
early explorations employed radiative-transfer simulations that did 
not include all the relevant physics, in particular lacked a rigor- 
ous treatment of non-LTE effec ts. In the context of SN 1987 A, 
lEastman & Kirshneil l ll989h and ISchmutz et al] ( 119901) were also 
unable to reproduce the He I lines at 1-2 d after explosion, and obvi- 
ously, non-thermal processes could not be the culprit. This He I dis- 
crepancy has now been resolved, although at a gr eat computational 
cost, with the use of a full non-LTE treatment jPessart & Hillieil 
llOlOl) . demonstrating the need for a fully-consistent solution of the 
radiative-transfer problem. More generally, it is interesting to in- 
vestigate whether non-LTE effects alone can lead to the produc- 
tion of He I lines in W-R-star explosions. Using the present simula- 



tions, we study this possibility at early-times when " Ni produced 
at depth has influenced the photospheric layers neither directly nor 
indirectly. We also include in this discussion those models that do 
not contain ^''Ni or other unstable nuclei. Our results thus apply 
irrespective of ^*Ni production and therefore yield stronger and 
un-compromised constraints on the progenitor star itself. 

The montage of spectra shown in Figs. |6}0 for 1.5 and 7.0 d 
post-explosion times suggest that He I lines are present, although 
sometimes quite weak, in all models apart from Smi60mf7p08zl. 
In model Bmi25mf7p3z0p2, the synthetic spectrum at 1.56 d is 
almost exclusively composed of He I lines, specifically at 3705, 
3819, 3888, 3964, 4026, 5875, 6678, 7065, 7281, 10830, 10914, 
and 1 1969 A (these are the He I lines with a Sobolev optical depth 
greater than 50 A in the range 3500-12000 A). When He I fines 
are rare and weak. He 1 5875 A is the most apparent in the op- 
tical range. Compared to what they were at 1.5 d, He I lines at 
7.0 d are of similar (and weak) strength in models BmiI8mf3p79zI, 
Bmil8mf4p41zl or Bmi25mf5p09zl, but of even lower strength in 
models Bmi25mf6p49zl and Bmi25mf7p3zOp2. This evolution is 
more pronounced for He 1 5875 A than for He 1 10830 A, which can 
remain relatively strong both in absorption and in emission up to 
15 d after explosion (Fig.|9). 

Overall, the presence of He I lines in our synthetic spec- 
tra suggests that He I is thermally and radiatively excited in 
the con'esponding models. The conditions are more favorable in 
progenitors with a sizeable and loosely-bound helium envelope 
(Fig. |4j, which tend to produce ejecta with a suitably larger in- 
ternal energy (for a given ejecta kinetic energy). Note that in 
some of our models, the feature at 5900 A is due to Na I D rather 
than He 1 5875 A (models Bmi25mf5p09zl, Bmi25mf6p49zl, and 
Bmi25mf7p3zOp2; Fig.^. 

In model Smi25mfl8p3zOp05, the only helium line present is 
He 1 10830 A, which appears as a weak ("detached") absorption at 
a Doppler velocity of ^^-26,000 km s~^, reminiscent of the line- 
profile structure found in time-dependent simu lations displaying 
an io nization freeze-out (see Figs. II and 12 of iDessart & Hillieil 
I2OO8I) . In this model, both ionisation freeze-out and chemical strat- 
ification are the cause of this feature. 
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Figure 9. Top: Same as Fig. [U but now for the Hel5875 A region (left; we also witness absorption from the 6678 A line to the red) and the Hel 10830 A 
region (right). The model depicted is Bmil8mf3p79zl. Middle: Same as top, but now for model Bmi25mf5p09zl. Bottom: Same as top, but now for model 
Bmi25mt7p3z0p2. In this low-metallicity model, Hel lines are present at early times, and persist even for two weeks for the neai'-IR line, even in the absence 
of non-thermal excitation by high-energy electrons. Note that the abscissa scale now extends to 50000 km s~^. 



In Ni-rich models, although we neglect high-energy elec- 
trons and treat all decay energy as a local heat source, He 1 lines 
can be present beyond the end of the post-breakout plateau. For ex- 
ample, He I lines are present up until peak brightness at ~30 d in 
model Bmil8mf3p79zl because at that time the helium mass frac- 
tion at the photosphere is still ~0.8 and the diffusing heat from de- 
cay enhances the gas temperature at the photosphere to 10000 K. 
In model Bmi25mf5p09zl, at 30 d after explosion, the helium mass 



fraction at the photosphere is down to ~0.5 and the heat from de- 
cay, shared by a greater mass, is also diffusing out more slowly 
(Fig. |2]l. This later point may be ciritical for the SN Ic classifica- 
tion: They may be genuinely helium-deficient, or may correspond 
to ejecta with a mixture of helium, carbon, and oxygen, and sep- 
arated from the Ni-rich layers by a sizeable mass buffer. In any 
case, we stress that, even as a pure heat source, decay energy dif- 
fusing out from greater depths may be sufficient to produce He I 
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Figure 10. Illustration of tlie mean helium ionisation level for models 
Bmil8mf3p79zl (top), Smi60mt7p08zl (middle), and Bmi25mf7p3zOp2 
(bottom). In each panel, the coloured dots refer to this ionisation level at 
the photosphere and its evolution with time. An interesting feature is the 
systematic partial ionisation of hydrogen and helium in models showing 
H I and/or He I lines at early times, even in the absence of ^^Ni/^^Co; the 
counterexample here is model Smi60mt7p08zl. [See text for discussion.] 



lines at the light-curve peak, provided the ratio of ejecta mass to 
^^Ni mass is not too large and the helium mass fraction at the pho- 
tosphere is high (as in this model Bmil8mf3p79zl). 

The notion that one needs He 1 10830 A to unambiguously 
assess the presence or absence of helium is partially supported 
by our simulations. At early times, whenever He 1 10830 A is 
present. He 1 5875 A is present too and should be identifiable. At 



~10d, however, He 1 5875 A tends to become very weak while 
He 1 10830 A may still be strong and more easily identifiable. 
He 1 10830 A is systematically broader and stronger, both in absorp- 
tion and emission, than other He I lines, and it suffers less from line 
blending. 

The disappearance of He I lines is directly connected to the 
decreasing helium mass fraction at the photosphere (Fig.|3), as it 
migrates towards the originally more tightly-bound inner regions 
of the ejecta. This phase also corresponds to the rapid luminosity 
fading that ends the post-breakout plateau (Fig.[T}. As for hydrogen 
discussed above, the helium ionisation is a critical component for 
the production of He I lines. In model Bmi25mf7p3z0p2, which 
shows the strongest He I lines at 1.5 d, helium is partially ionised 
above the photosphere, while in other simulations with weaker He I 
lines, it is very nearly neutral (Fig. 110). 

There is something quite spectacular here though. While we 
find that hydrogen present at a few percent by mass is sufficient to 
produce a very strong Ha line, helium present at ~90% by mass 
may produce feeble He I lines, as in simulations Bmil8mf3p79zl, 
Bmil8mf4p41zl or Bmi25mf5p09zl. It is not clear at what thresh- 
old helium mass fraction He I lines would disappear; in our re- 
stricted set of simulations. He I lines are invisible for helium mass 
fractions of <50%. In these models, the lines we do predict stem 
from species that have a mass fraction at the per cent level or less. 
This is a dramatic illustration of the direct effects of ionisation and 
excitation on the emergent radiation, whi ch can easily compensat e 
for even large variations in abundances dOessart & Hillieij|2008h . 
Accounting for such "dark matter" is a challenge. Of course, the 
presence (and dominance) of helium does condition the radiative- 
transfer solution (for example, helium ionisation controls the loca- 
tion of the photosphere). It also conditions the expansion rate since 
helium is associated with a sizeable mass, but its total mass cannot 
be derived easily in a direct manner. 

All our binary-star evolution models are characterised by very 
high helium-surface mass fractions of >85% and their associated 
ejecta produce spectra that contain He I lines (potentially difficult 
to identify in the optical but not with He 1 10830 A). The identifi- 
cation of He I lines prior to the SN re-brightening phase is a criti- 
cal signature of a very high helium content, potentially favouring a 
binary-star evolution channel. Our results also imply that a signif- 
icant amount of helium may be present in the outer regions of the 
ejecta without any clear radiative signature. Recall that we are ig- 
noring non-thermal processes, which could help produce He I lines 
when we do not predict them. However, in those cases where we do 
predict them, such processes would only make He I lines stronger. 



4.3 CNO elements 

With more erosion of the progenitor star and/or helium core, deeper 
shells are revealed and one expects a progression from a nitrogen- 
rich (and potentially hydrogen-rich) to a carbon-rich, and then to 
an oxygen-rich W-R star. In the same sequence the helium mass 
fraction should steadily decrease. Our set of progenitor models re- 
flects this dichotomy. The binary-star models correspond to WN 
stars since they show enhancements of helium and nitrogen (with 
possible traces of hydrogen), and depletion of carbon and oxygen, 
as expected for the CNO cycle. The single-star models correspond 
to WC or WO stars since they show depletion of helium and nitro- 
gen (with no trace of hydrogen) and enhancements of carbon and 
oxygen, as expected from helium-core burning. Hence, variations 
in hefium mass fraction are accompanied by coiTesponding varia- 
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tions in CNO abundances. It is tiius important to search for radiative 
signatures of sucii CNO elements. 

As expected, those (single-star) models that do not show 
any He I line for the first 10 d show a wealth of lines from C I 
and O I in the red part of the optical range, while only very few 
and weak lines are associated with nitrogen (i.e., Nl; Figs. [6}- 
[TJ. There are numerous C I lines, but the stronger ones are at 
8335, 9405, 9658, and 10691 A. For Ol, the stronger lines are at 
7774, 7987-7995, 8446, 9261, and 11287 A. Some of these over- 
lap with He 1 10830 A, but He 1 5875 A should then be searched 
for in early-time observations to help resolve the ambiguity. There 
are also lines from once-ionised species such as C II, for exam- 
ple in model Bmi25mfl8p3zOp05 at 1.51 d, with the strongest fea- 
tures at 5890 and 6580 A (doublet lines; not shown here). All these 
C/O lines have similar strength in models Smi60mf7p08zl and 
Smi25mfl8p3zOp05. The appearance of such C/O lines arise from 
the large C/O mass fraction at the photosphere in both models, each 
element having a mass fraction of 0.2-0.5 (Table Such single- 
star progenitors would be classified WC stars and the lack of helium 
lines in their optical spectra at 1-10 d after explosion, caused by an 
ejecta helium mass fraction <50%, would make the coixesponding 
explosion look like a Type Ic SN. 

In contrast, models for which we predict He I lines do not 
show strong C/O lines. This occurs because of the lower photo- 
spheric C/O mass fraction. The nitrogen mass fraction is also low, 
but not as low, and so we predict some N I lines in the red part 
of the spectrum. In simulation Bmil8mf3p79zl at 7.0 d, there are 
numerous Cl, N I, and O I lines of weak/moderate strength. The 
main N I line is at 8680 A, and as we shall see further below, and 
from what we discussed above, all these lines/species contributions 
overlap and challenge an accurate identification. 

4.4 Intermediate-mass elements 

At the early times, the photospheric composition reflects only the 
results of (steady) core burning, causing modest enhancements of 
sodium, neon or magnesium abundances compared to solar. In con- 
trast, titanium or silicon are at the environmental abundance, just 
like iron. Despite their generally low abundance, a number of IMEs 
have associated line features, in particular Na, Mg, Si, and Ti. This 
is analogous to what is observed in Type II SNe where species at the 
envrionmental metallicity cause significant blanketing (e.g.. Till), 
or produce obvious features associated with resonance lines (e.g., 
Nal or Call). If we evolved these ejecta to later times, lines from 
such IME would be favoured by the increasing mass fraction of 
such elements at the photosphere. 

NalD is weak or absent at 1.5 d, but strengthens as the pho- 
tosphere cools. At 7.0 d, it appears as a strong line in all our mod- 
els, whether at solar or sub-solar metallicity, and thus even for an 
Na mass fraction of 10^''. In models showing helium lines, we find 
that it overlaps with He 1 5875 A and contributes significantly to the 
feature. 

Our single-star models, which show numerous and strong C I 
and O I lines, also show strong Mg II lines, primarily at 4481, 7877- 
7896, 8213-8234, 9218-9244, 9632, and 10915 A. 

Irrespective of abundance issues, Ca always shows strong 
Ca II lines as soon as it becomes once ionised. In all our simula- 
tions, we predict a strong Ca II H&K as well as a strong Ca II triplet 
at 8500 A. This holds even in low-metallicity models once they pos- 
sess the required ionisation. 

We also predict a few lines of Si, in particular Si II at 6347- 
6371 A. In the model Bmil8mf3p79zl, it contributes a line that 



overlaps quite closely with Ha, although it is weaker than Ha. In 
model Bmi25mf5p09zl, this Sill doublet appears as a strong P- 
Cygni profile at 1.5 d, and could fuel some controversy about the 
presence of hydrogen in the ejecta (present with a mass fraction 
<10~^). As we demonstrated before, early-time observations can 
help lift this ambiguity, since even very small amounts of hydrogen 
would be sufficient to produce a strong, unambiguous. Ha line. 
Models in which Ha is strong at 1.5 d also show strong Ha at 7 d. 
In contrast, the Sill doublet in Bmi25mf5p09zl is strong at 1.5 d 
but has vanished by 7.0 d after explosion. 

4.5 Iron-group elements 

The main iron-group element that contributes strong line- 
blanketing is iron (the IME titanium also causes significant blan- 
keting, primarily in the B band). Due to the low photospheric tem- 
peratures, Fe II is the dominant ion and causes numerous absorption 
(and some emission) lines throughout the optical range. This occurs 
in all solar-metallicity models presented [i.e., a large iron abun- 
dance is not required to produce significant iron line-blanketing 
jHarkness et al.ll9^7h l. However, a small iron abundance, as in the 
low-metallicity model Smi25mfl8p3zOp05, considerably reduces 
the strength of line blanketing. In that model at 6.93 d, the spectrum 
shows only very weak Fe II lines, in contrast with the dominating 
C I/O I/Na I/Mg Il/Ca II lines. 

This is a very interesting result which suggests that quanti- 
tative spectroscopic analyses of SN spectra may, when they have 
become accurate enough, be a tool for the determination of the SN 
environmental metallicity. This will require early-time observations 
since it is only at such epochs that we are confident the photosphere 
is only a probe of the progenitor surface layers, uncormpted by ele- 
ments synthesised explosively at depth. However, it may be limited 
to events in which external disturbances are negligible (i.e., cases 
where interaction with the pre-SN wind or with a companion star 
is weak). 



5 COMPARISON WITH OBSERVATIONS AND 
DISCUSSION 

In the preceding section, we have presented a number of spectro- 
scopic signatures for ejecta stemming from core-collapse explo- 
sions of W-R progenitor stars. Let us now discuss the main results 
and confront them to observations. 



5.1 Post-breakout plateau luminosity 

A generic feature of all synthetic bolometric light curves pre- 
sented in this work is the presence of a post-breakout plateau 
brightness that typically lasts for ~IOd (Section [3j see also 
[Ensman & Wooslev 1988). The origin of this plateau is analo- 
gous to what is seen in Type II SNe, but differs in brightness be- 
cause of the small progenitor radius, the smaller ejecta mass, and 
progenitor composition. In Fig. [TT] we illustrate the diversity in 
SN light curves with different ejecta models all having 1-I.2B 
kinetic energy. We consider three cases — progenitor stars hav- 
ing a small size (~12.3R0; model Bmil8mf4p41zl, this work), 
a modest size (~50Rp) ; model "lmI8a7Ad" for SN 1987A from 
iDessart & Hilliej|2010l but now evol ved to nebular times), a nd a 
large size (~800Ro; model sI5el2 of lDessart & Hilliej20Ilh . All 
show a post-breakout plateau which tends to be brighter for larger 
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Figure 11. Comparison of the bolometric-luminosity evolution of SN lb 
model Bmil8mf3p79zl (0.184 M0 of ^'^Ni), SN Il-pec mo del "ImlBAf 
of SN 1987A (0.084 Mq of ^SNi; lOessart & HUUeill2010l) . and SN Il-P 
model sl5el2 (0.086 of ^"^Ni: IPessart & HiUieill201ll) . These three 
models have a similar ejecta kinetic energy of 1-1. 2B, and full 7-ray trap- 
ping is assumed (it does not hold for model Bmil8mf3p79zl past >70d; 
see Fig. 113) . We also overplot the power associated with the decay of (ini- 
tially) 0.05 Mq of ^''Ni (olive curve). All models show an early-time lumi- 
nosity plateau whose magnitude is function of the progenitor radius and the 
H/He abundance ratio, and a late-time luminosity behaviour reflecting the 
initial amount of ^®Ni. The distinct evolution is primaiily conditioned by 
the progenitor-envelope properties, i.e., in paiticular whether the progenitor 
radius is ~1OR0 (Bmil8mfip79zl), ~5OR0 (lml8a7Ad), or ~8OOR0 
(sl5el2). 
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Figure 12. Ejecta-temperature distribution versus radius at ~10d af- 
ter explosion for the SN Ic model Smi60m t7pO8zl, SN lb mode l 
Bmil8mf4p41zl, SN ll-pec m odel lml8a7Ad ^Dessart & HillieJ I2OIOI) , 
and SN 11-P model sl5el2iso )Dessart & Hilliej|201ll) . We show CMF- 
GEN results (solid lines; the photosphere location is shown as a dot), as 
well as results from a separate code that assumes no diffusion of heat and 
pure adiabatic cooling (broken line; our "1-zone model"). The contrast in 
brightness between models (Figures [T] and II U is reflected by the larger 
radii/temperatures of the corresponding ejecta and epoch, and nearly ex- 
clusively reflects modulations in (adiabatic) expansion losses. 



progenitor objects, and in the range 4x 10^ to 6x 10* L©. Its du- 
ration is conditioned by the progenitor-envelope binding energy, 
ejecta mass, composition, as well as the mass and distribution of 
^®Ni. Here, it varies from a few days in model Bmil8mf4p41zl, 
to about 20 d in mo del "lml8a7Ad" (co mpatible with the observa- 
tions of SN 1987A. |Phillips et al.lll988l) . and about 100 d in the SN 
II-P model sl5el2. 

This light-curve diversity is also directly connected to the evo- 
lution of the ejecta internal energy with time, which is drastically 
different between such ejecta. To emphasize the importance of adi- 
abatic cooling through expansion, we have evolved a number of 
SN ejecta with a separate program that assumes no diffusion of 
heat, pure adiabatic cooling, and treats radioactive-decay energy 
as a pure local heating source. In this respect, it evolves each ejecta 
mass-shell individually (i.e., one-zone model), updating its energy 
according to Ae = -pA(l/p) + (Jedocay, where e, 5edccay, P, 
and p are the specific energy, the decay energy, the pressure and 
the density associated with that mass shell. We solve this equation 
using the Newton-Raphson technique and employ an equation o f 
state, developed as part of the work shown in lOessart et al. I dlOlObh , 
which is function of density, temperature, and composition. 

As shown in Fig. [T2] for a post-explosion time of ~10d, 
the inner-ejecta temperature obtained with CMFGEN agrees very 
well with the one computed assuming no diffusion and pure adi- 
abatic cooling (this also serves as an independent and additional 
check on CMFGEN predictions). Because of the huge range in 
progenitor radii, these inner-ejecta temperatures vary by up to a 
factor of ten between the ^^Ni-deficient originally-compact SN 
Ic model Smi60 mf7p08zl and the origin ally-extended SN II-P 
model sl5eI2iso jPessart & HillieJ I2OI ll) . Closer to the photo- 
sphere (shown as a dot), the temperature is influenced by non-LTE 
effects (sl5el2iso), radiative cooling (lml8a7Ad), and the heat 
wave generated by decay heating at depth (Bmil8mf4p4Izl). The 
contrast in brightness between models (shown in Figures [Tlandlllb 
is reflected by the larger radii/temperatures of the corresponding 
model and epoch, and nearly exclusively reflects modulations in 
(adiabatic) expansion losses. 

The post-breakout plateau is generally only observed in SNe 
II, and in particular those of the plateau type because it is bright 
and lasts for a long time. In contrast, SNe Ib/c are generally dis- 
covered during the re-brightening phase, around peak light, or even 
later, so that any plateau phase taking place earlier is missed, and 
the shock-breakout time is poorly constrained. There are a few ex- 
ceptions to this. For example, the fading from the breakout phase 
was observed in SN 1 993 J, which transit ioned after a few day s 
to a re-brightening, but by only ~1 mag l lRichmondetal]|l994h . 
SN 2005bf was a peculiar SN lb in that it showed a double-peak 
bolometric light curv e. No plateau is visible prior to the first peak 
dFolatelli et al.l2006h and the explosion time is not known, although 
this SN was not discovered earlylf] In cont rast, SNe associa ted with 
7-ray-b ursts (GRBs), e.g., SN I998bw jPatat et al.ll20oTh or SN 
20I0bh dChomock et al.l2010l) . have a well defined explosion time, 
but in these, the object seems to be always brightening at >I d after 
the GRB signal . SN 2008D was caugh t as the shock broke out of the 
progenitor star ([Soderberg et al. l2008l:IChevalier& Franssonl2008l ; 



iModiaz et alj2009l) . By I d after breakout, the luminosity plateaus, 



^ The first peak shows a rise time and a maximum luminosity that is in fact 
very compai'able to standard SNe lb. The "anomaly" with SN 2005bf is the 
bright second peak, which does not seem to stem from decay energy but 
seems more consistent with the birth of a magnetar l lMaeda et alj2007h . 
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but already by 5 d the SN re-brightens. Interestingly, it brightens by 
merely 1 mag to reach a peak before fading again. 

The post-breakout plateau that we discuss here seems to have 
been seen in SN 2008D, but lasts for a few days only and has 
a a brightness that is merely 1 mag below the peak. This small 
post-breakout-plateau/peak brightness contrast is in part caused by 
the faint light-curve p eak, supporting a smaller-than-averae e ^"Ni 
mass of ~0.07Mq ffanaka et al.ll2009l ; iDrout et al.l 12010) . com- 
pared to ~0.2 M0 in our models. The short post-breakout plateau 
suggests that the ejecta mass is small and/or mixing qu ite efficient. 
The huge plateau brightn ess of ~1.5x10**Lq tSoderberg et al.l 
l2008l : iModiaz et aU 120091) is, however, perplexing. This plateau 
brightness is thrice that obtained in our simulations. It cannot be 
due to a lar ge progenitor radius since analysis of the early X-ray 
light curve jSoderberg et al.ll2008l: iModiaz et aLll2009t) supports a 
progenitor radius of i?* <J 1 R0. Efficient outward mixing of '''^Ni 
would enhance the brightness, but it would likely induce a bright- 
ening rather than a plateau, and is therefore not the likely explana- 
tion. Some light contamination could come from the galaxy host, 
although it appears unlikely here since SN 2008D was followed up 
to very late times when the SN was visually much fainter than im- 
mediately after breakout. The contamination could come from the 
SN light itself, but emitted earlier and scattered by the surrounding 
CSM or pre-SN mass loss. The lack of narrow lines from photo- 
ionised and/or shocked CSM gas in SN 2008D early-time spectra 
does n ot strongly support this. Furthermore, Chevalier & Franssor3 
( l2008l) argued for a low mass rate for the SN 2008D progenitor W- 
R star. Following upon the argument that the SN 2008D progenitor 
star may be in a binary system, the large early post-breakout lumi- 
nosity could in part arise from th e collision of the SN ejecta with 
the companion star ( lKasen|[201(t) . The same argument could also 
explain the early light curve of SN 1993 J, which in current mod- 
els requires a low-mass H-poor and He-rich p rogenitor with a sur- 
prisingly large radius of ^630 R© i jBlinniko'v et al. 1998). Early- 
time light curves, starting immediately after breakout, may there- 
fore provide an important clue to the single/binary status of the 
progenitor star, and complement independent constraints set by the 
outer-ejecta composition (Sectionl?^. 

Compared to the progenitor W-R stars, which emit the bulk 
of the radiation in the UV, the p ost-breakout plateau should be 
> 10 m ag brighter in the \/-band ( I van der Huch3l200ll : ICrowthed 
l2007l : this estimate accounts for the large bolometric correction of 
W-R stars). This enormous visual brightening might still be difficult 
to detect if the SN is located in a crowded region or if it falls on the 
galaxy light. For SNe located at large galactocentric radii, it should 
however be possible, provided the search is deep and performed on 
a daily cadence. 

5.2 Peak luminosity and width 

In Nature, '''^Ni mixing within SN Ilb/Ib/Ic ejecta will be more effi- 
cient than in our set of ^''Ni-rich models, which assume no mixing 
at all. Our neglect of non-thermal processes will affect the gas ioni- 
sation and compromise our synthetic colors at peak, which are thus 
not discussed in this comparison to observations. We thus discuss 
instead the bolometric properties. 

Our ^'^Ni-rich models Bmil8mf3p79zl and Bmil8mf4p41zl 
show a 20-d rise time from the post-breakout plateau to the peak, 
whose bolometric magnitude is -17.89 and -17.64 and bolomet- 
ric luminosity of 1.14x10® and 9.0xI0*Lq, respectively. These 
trends reflect the larger ^®Ni mass in the former model (0. 1 84 com- 
pared to 0. 170 Mq ) combined with the larger ejecta mass in the lat- 



ter (2.39 compared to 2.91 Mq). Importantly, ' Ni is present out to 
the same ejecta velocity in both. 

In contrast, the bolometric light curve for model 
Bmi25mf5p09zl has a longer rise time of 40 d from the end 
of the post-breakout plateau, with a peak bolometric magnitude 
(luminosity) of -17.49 (7.9xlO®L0). Here, despite the larger 
^^Ni mass of 0.24 M©, the peak is fainter than for the other 
two ^^Ni-rich models because of the larger ejecta mass and the 
significantly deeper location of ^"^Ni, i.e. below 1250km s~^ 
compared to below ^2500 km s~^ . 

Furthermore, a larger ejecta mass and a deeper location 
of ^®Ni broadens the peak width significantly. Fifteen days af- 
ter the light curve peak, we obtain a magnitude fading of 0.82, 
0.57, and 0.15 for models BmiI8mf3p79zI, Bmil8mf4p41zl, and 
Bmi25mf5p09zl. Modulations in ejecta mass and ^''Ni mixing can 
be interchanged to yield similar results. In a separate set of simula- 
tions for these three progenitors but characterised by a ^''Ni extent 
out to 3500 km s^^, all three models yielded similar light curve 
widths. In SNe Ibc light-curve modelling, mixing is usually ad- 
justed freely to obtain a good match but we see here that there is a 
fundamental degeneracy in the problem. This requires an alternate 

constraint on mixing, which can only come fr om spectra. 

Overall, our results agree with those of lEnsman & WooslevI 

( Il988l) when adopting similar ejecta characteristics. Our main char- 
acteristics for models Bmil8mf3p79zl and Bmil8mf4p41zl are in 
good agr eement with t he ob servations and estimates recently pub- 
lished bv lDrout etai] ( l2OI0h and b ased on a large sample of SN e 
Ibc multi-color light curves (see also lRichardson et al.l2003.l2006h . 

5.3 Post-peak luminosity fading 

The SN post-peak/nebular luminosity is understood to be powered 
by radioactive decay, and primarily that of ''^Co at the times of in- 
terest here (we consider explosions that synthesise primarily '''^Ni). 
In SNe II, the large ejecta mass of typically 10-15 M© ensures the 
full trapping of 7-rays for a few hundred days. However, in Type I 
SNe, the lower ejecta mass is less efficient at trapping 7-rays, which 
can escape e arlier (for a similar phenome non in the context of SNe 
la, see, e.g., iHoeflich & Khokhlo\i[l993) . How much earlier, and 
how efficiently, are a very important question, since it connects 
to the ejecta mass and location of ^^Ni. The post-peak/nebular- 
phase luminosity decline rate represents an important diagnostic 
for the explosion and the progenitor star. Because 7-rays are in- 
sensitive to the gas ionisation state and composition (provided only 
elements heavier than helium are present), 7-ray transport is con- 
siderably simpler than optical-photon transport and less subject to 
uncertainty. 

The majority of SNe Ilb/Ib/Ic show post-peak luminosity de- 
cline rates that are st eeper than the expected r ate of O.OI mag/d for 
full 7 -ray trapping dRichardson et al.l l2006f) . lEnsman & WooslevI 
( Il988l) invoked the effect of clumping. Departures from sphericity 
or efficient in trinsic mixing may als o propel ^^Ni further out (see, 
for example, iTominaga et al.l l2005l : iFolatelli et al.l |2006|) . which 
would then increase 7-ray escape during the nebular phase. Un- 
fortunately, the emblematic SN 1987 A, for which there is strong 
evidence of an asymmetric explosion and for which early 7-ray es- 
cape is inferred, shows a nebular-luminosity fading rate that ex- 
actly matches the O.OI mag/d decline rate that implies full trap- 
ping jAmett et al.l 19891 ). As alluded to above, a systematically low 
ejecta mass combined with moderate mixing may alone explain the 
generically fast decline rate of SNe Ilb/Ib/Ic. 

In the CMFGEN simulations we have presented so far, ^''Ni de- 
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Figure 13. Evolution of the bolometiic luminosity of our non-LTE time-dependent radiative-transfer simulations (Lbob sohd lines), which assume full- 
trapping of 7-ray photons, and based on the ^^Ni-iich models Bmil8mf3p79zl(black), Bmil8mf4p41zl (blue), and Bmi25mf5p09zl (red). The dots con'e- 
spond to the actual epochs computed. If we had instead used in CMFGEN the 7-ray-energy deposition function output by our 7-ray-Monte-Carlo transport code, 
the bolometric luminosity would have turned over at ~70d in models Bmil8mf3p79zl and Bmi25mf4p41zl and fallen-off more steeply than the 0.01 mag/d 
rate for full-trapping of Co-decay photons (it would have followed the dashed curve representing the energy effectively deposited per unit time in the ejecta 
'^edcp.dccay/t'i)- The dot-dash line corresponds to the total energy released, in the form of 7-rays and positrons, by radioactive decay in the Monte Carlo 
simulations (it is offset by <10% from the CMFGEN curves at very late times because of the numerical diffusion inherent to the necessary regridding at each 
time step in CMFGEN). The smaller the mass overlying the "® Ni-production site, the steeper is the luminosity-decline rate. A low ejecta mass may thus be one 
way of explaining naturally the rapidly fading light curves of most SNe Ilb/Ib/Ic, without invoking the intervention of clumping, extreme mixing, or jets. 



cay energy, when included, was fully deposited at the site of release 
in the ejecta. With this full-trapping assumption, our corresponding 
models fade at a 0.01 mag/d at nebular times. To check how the 
light curve may be affected by this approximation, we have devel- 
oped a 7-ray-transport Monte Carlo code that computes the fraction 
of 7-rays that escape, and of those that are trapped, how they dis- 
tribute their energy within the ejecta (Hillier & Dessart, in prepa- 
ration). In Fig. [T3] we show results based on our *'^Ni-rich ejecta 
models Bmil8mf3p79zl, Bmil8mf4p41zl, and Bmi25mf5p09zl 
(and assuming spherical symmetry). We find that the smaller the 
ejecta mass (for the same ejecta kinetic energy, here of 1.2B), the 
earlier the fading rate departs from the rate for full-trapping and the 
steeper the luminosity fades subsequently. 

In models Bmil8mf3p79zl and Bmi25mf4p41zl, the post- 
peak luminosity curve smoothly joins the steep nebular decline, be- 
cause the 7-rays start escaping immediately at the fall off from the 
peak0 The subsequent fading is very steep, and compatible with 



^ Our 7-ray-trans port Monte Carlo simu lations for SN 1987A based on 
model "ImlSAf loessart & Hilliej|2010h predict the onset of escape for 
the bulk of 7-rays at >500d after explosion. The st eeper decline of th e 
SN 1987A bolometric luminosity occurs at that time jAmett et alll 19891) . 
but un fortunately this is co ncomitant with dust formation in the SN 1987A 
ejecta )Colganetai]|l994) . which may be responsible for part of the in- 
ferred attenuation. 



the decline rate of SNe Ib/Ic jRichardson et al.ll2006l) . However, in 
model Bmi25mf5p09zl, full-trapping persists at the onset of the 
nebular phase and it is only later that 7-ray escape occurs and be- 
comes noticeable. There is thus a point of inflection at 130d in the 
bolometric luminosity when full-trapping ceases to hold. Observa- 
tions do not seem to support such late trapping, but instead suggest 
at least some moderate level of mixing. In model Bmi25mf5p09zl, 
adopting some mixing of '■'■'^Ni out to 3000 km s"^ (instead of 
1250km s"^) moves the onset of 7-ray leakage to 70 d after ex- 
plosion. 

In the observa t ions of SNe 19 93 J jUtrobinl 1 1994 
Wooslevetal] Il994j; lYoung et alj Il995l). 19941 jYoung et'all 



19951 : 



Iwamoto et al 



1994h . or 1998bw jWooslev et al.ll 19991) . the 



post-breakout re-brightening occurred very soon after breakout, 
the peak was reached very early, the rise and fall from the peak 
occurred over a short time, and the decline rate was faster than 
for full-trapping. All these independently support the notion, 
which is not new, that SNe Ilb/Ib/Ic are associated with low-mass 
ejecta, and quite generally low-mass progenitor stars. However, 
it may be that rather than invoking extreme clumping or mixing 
configurations, one may be able to accommodate such observations 
with moderate mixing and low-mass ejecta exclusively. In the 
context of hypemovae, associated or not with a 7-ray signal, 
the larger ejecta expansion rate reduces the diffusion time and 
enhances the probability of 7-ray escape, so that a larger mass may 
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still be compatible with the observations. In the future, to address 
this issue, we will model simultaneously, in the fashion described 
here, the spectra and light curves of these rare events. 

5.4 Early-time spectra 

Our early-time synthetic SEDs look like that of a cool thermal emit- 
ter, with a peak distribution around 5000 A, and affected by nu- 
merous weak/strong lines from H/He/C/N/0 and IME as well as 
forests of lines from blanketing species and in particular iron. The 
main results from our work are that hydrogen and/or helium lines 
are predicted even in the absence of non-thermal processes, or heat 
diffusion, caused by radioactive decay. Hi Balmer lines are seen 
for a hydrogen mass fraction as low as 0.01 in the corresponding 
ejecta layers, while He I lines seem guaranteed only if the helium 
mass fraction is close to 100%. During the rise to peak, decay heat- 
ing may be strong enough to produce He I lines, which may persist 
until the peak of the light curve in spite of the much lower helium 
mass fraction at the photosphere at that time. Such early-time ob- 
servations (spectra and light curves) are critical since during the 
post-breakout plateau, up to <50% of the ejecta mass may have 
gone through the photosphere, and thus potentially be unaccounted 
for. 

There is a scarcity of SN Ilb/Ib/Ic spectroscopic observations 
at early times, prior to the post-breakout re-brightening. Indeed, 
the first spectra are typically taken during the rise to peak, which 
is a phase influenced indirectly by radioactive decay (an outward- 
diffusing heat wave). The only exceptions with early-time observa- 
tions are SNe 1993 J (lib) and 2008D (lb). Numerous investigations 
have studied the light curve or the spectra, but none have performed 
the full non-LTE time-dependent radiative-transfer simulation of 
the ejecta, with the approach presented here. It is critical to do both 
at the same time, to test that the ''''Ni mixing adopted to fit the 
light curve is also compatible with the excitation/ionisation seen in 
spectra. This has never been done. 

Our models were not tailored to match any specific obser- 
vation so we cannot be quantitative. Qualitatively speaking, our 
model Bmil8mf3p79zl has some analogy with the early obser- 
vations of SN 1993J (e.g., the presence of a strong Ha line). We 
note that the ea rly- time SN 1993 J spe ctra are nearly featureless 
and rather blue dMatheson et al.ll2000l) . As discussed above, this 
has been associated in former studies with the exp losion of an ex- 
tende d low-mass progenitor in a binary system dSlinnikov et al.l 
Il998h . Alternatively, the large early-time flux could in a large part 
stem from the collision of the SN ejec ta with the identifi ed com- 
panion of SN 1993J's progenitor star dMaund et al.ll2004 . which 
would call for a significant revision of current models of the event. 
In model Bmil8mf3p79zl, the pre-SN star has a radius of 10 Rq 
and the early-time spectra are consequently rather red (the photo- 
sphere is cool and recombined), not blue, and line blanketing is 
strong. Irrespective of the mismatches with observations, we do 
concur with previous investigations that low-hydrogen mass frac- 
tions and low-hydrogen integral mass at the progenitor surface are 
sufficient to produce a strong (and broad) Ha line in t he corre- 
sponding early-time SN lib spectra, as observed in 93J ( UtrobinI 
_1 994 ; 'Wooslev et al. [ll994j : lYoung et al.lll995l : iBaron et al.lll995l : 
Blinnikov et al. 199: 



None of our models match the observations of SN 2008D dur- 
ing the first week after explosion. Given our results for hydrogen, 
it is clear that 08D is hydrogen-deficient. However, at early times, 
the observed spectra are nearly featureless, as for SN 1993J, while 
our simulations predict He I or Fe II lines for such a SN lb event. 



One striking property of SN 2008D is that the early post-breakout 
plateau is about three times as bright as that of SN 1987 A and 
only < 1 mag fainter than the light-curve peak. SN 2008D is much 
brighter at those early times than expected for a rather compact 
W-R-star progenitor Together with the nearly featureless spectra, 
these properties suggest that the SN light is indeed "corrupted". As 
discussed above, an attractive solution to this is emission from the 
collision of the SN ejecta with a companion star ( .Kasen.2010) . 

At the time of re-brightening in model BmiI8mf3p79zI, the 
heat wave causes the appearance of a kink in the Ha absorption. If 
stronger, the corresponding change in ionisation could have led to 
a reversal of the Doppler velocity of maximum absorption in Ha. 
In model Bmi25mf5p09zl, the effects of the heat wave is very pro- 
nounced and causes the photosphere to move out mass, or velocity, 
space. However, it occurs too late to leave any imprint onto He I 
lines. The velocity reversal in the He 1 5875 A line of SN 2005bf 
jModiaz et al.ll2009l) may stem from a similar circumstance, but in 
this case at a time when the photosphere is close to helium-rich lay- 
ers. T he origin of this uni que observation may be the birth of a mag- 
netar jMaeda et al1l2007l) . In general, such reversals/kinks/notches 
are not observed because SNe Ilb/Ib/Ic are discovered too late (dur- 
ing the rise to the peak), at a time when the photosphere is al- 
ready in the vicinity of the core. It is a dramatic change in heat- 
ing/excitation/ionisation conditions that permits this reversal, such 
as when the photosphere starts feeling the energy from decay or, in 
exceptional circumstances, from the magnetar radiation. 



5.5 ^ Ni-deflcient W-R-star explosions 

Our ^"^Ni-deficient models cannot be compared to observations 
since all SNe Ilb/Ib/Ic are discovered by the virtue of being ^®Ni 
rich. However, this could stem from the current bias for the discov- 
ery of brighter objects in bright galaxies, making objects devoid 
of a post-breakout re-brightening impossible to detect. Future deep 
and un-targeted surveys may in fact discover ''^Ni-deficient core- 
collapse SN explosions, and if they do so, it will be either through 
the breakout signal or the post-breakout plateau. As discussed in 
this paper, the radiative signatures during that phase are conditioned 
by the properties of the progenitor outer envelope. If detected, they 
would provide important information on the progenitor and per- 
haps clues on what distinguishes progenitors that explode as bright 
or faint supemovae, despite having a standard ejecta-kinetic energy 
of IB. 



6 CONCLUSIONS 

We have presented non-LTE time-dependent radiative-transfer sim- 
ulations of SN ejecta, with and without ^^Ni-decay products, re- 
sulting from the core-collapse explosions of single and binary 
W-R stars evolved at solar and sub-solar metallicity. Our ap- 
proach allows the simultaneous computation of the spectra and 
light curves, and is comp arable to our ea rlier studies of SN 1987 A 
jPessart&HiiiieJ 120101) and SNe II-P i Dessart & Hime3l201ll) . 
We pay particular attention to binary-star models for the produc- 
ti on of SNe Ilb/Ib/Ic , using the physically-consistent calculations 
of lYoon et al. U2O10h . This work contains about 400 separate cal- 
culations performed in separate time sequences, each time taking 
typically 1-2 days of computing, from which we infer the follow- 
ing: 

(i) All our SNe go through a ~10-d-long post-breakout plateau 
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with a luminosity of l-5xlO^L0, which is optically ~ 10 mag 
brighter than its progenitor In models endowed with 0.17-0.24 Mq 
of ^^Ni initially, this plateau is followed by a 20-30 d re-brightening 
phase up to a peak luminosity of 8 to 10 x 10* L©. Assuming full 
7-ray trapping, the nebular-luminosity decline rate of such ejecta 
agrees with the expected rate for ''''Co decay. However, 7-ray es- 
cape in our low-mass ejecta with moderate ^''Ni mixing becomes 
effective at ~70 d after explosion and lead to a much faster fading 
rate of ~0.02 mag/d. 

(ii) The absence of a sizeable post-breakout plateau, the narrow 
peak width, and the fast nebular-luminosity decline rate of most 
observed SNe Ilb/Ib/Ic light curves may be explained by a revi- 
sion downward of their associated ejecta mass (as well as of the 
ejecta kinetic energy to retain the same ejecta kinetic energy per 
unit mass)[3| This may be more physically justifiable than invoking 
strong clumping, efficient mixing to the progenitor surface, or ex- 
treme asymmetry of the explosion in the form of a "jet". Observing 
the post-breakout plateau and the post-peak luminosity decline rate 
of all SNe Ilb/Ib/Ic is key to further constrain this issue. 

(iii) We find that by the time of re-brightening, up to half the 
ejecta mass has gone through the photosphere. Studies on SNe 
Ilb/Ib/Ic based on peak and post-peak observations exclusively are 
thus subject to major shortcomings, since they miss important in- 
formation on the outer-ejecta properties, and in particular the H/He 
composition and the total ejecta mass. 

(iv) We present results for our W-R explosions with standard 
ejecta kinetic-energy but no explosively synthesised ^''Ni. These 
models go through a similar post-breakout plateau but irrevocably 
fade away as the photosphere recedes to helium-deficient inner re- 
gions, at about 10 d after explosion. 

(v) Hydrogen present with a small mass fraction of ~0.01 and a 
small cumulative mass of ^0.001 Mq gives rise to H I Balmer lines 
and in particular strong Ha in our WN-star progenitor models over 
the period 1-lOd. As latter times only Ha remains. However, it is 
weak and broad which, combined with possible overlap by C II and 
Si II lines, makes its identification difficult. Observed spectroscop- 
ically a few days after shock breakout, these events would even- 
tually be classified a s Type lib SNe (and correspond to the SNe 
cllb of lChevalier & S oderberg 2010), and probably as Type lb SNe 
otherwise. To address this issue, it is critical to obtain early-time 
spectra, preferably before the time of re-brightening (which may 
be caused by the recession of the photosphere to deeper, possibly 
H-deficient, ejecta layers). Using spectra taken at the peak of the 
light curve can hardly help resolve this problem since the photo- 
sphere may be a solar mass below the progenitor surface at that 
time. 

(vi) Simulations based on binary-star W-R (i.e. WN) progeni- 
tor models that have an helium mass fraction of ^95% in the outer 
~ I Mq of the ejecta show He I lines throughout their post-breakout 
plateau phase, irrespective of the absence/presence o/^^Ni. Past the 
end of the plateau, some models still show He I lines, even though 
our decay energy is treated as a pure local heat/thermal source. 
Non-thermal excitation/ionisation is however key to strengthen He I 
lines further, as observed. In general, the He 1 10830 A is stronger, 
broader, and longer lived than other He I lines like He 1 5875 A. 



^ A lower mass eje cta may help resolve the discrepant fitting of the SN 
1998bw light curve l lWooslev et al.|[l999l) , although it might then compro- 
mise our current understanding of the progenitors of long-soft GRBs. In the 
future, we will specifically study the case of high-energy W-R stai' explo- 
sions yielding ~10B-kinetic-energy ejecta to address this issue. 



(vii) Simulations based on single-star W-R (i.e. WCAVO) pro- 
genitor models that have an helium mass fraction of <50% in the 
outer ^ 1 Mq of the ejecta do not show He I lines during their post- 
breakout plateau phase, but instead lines of C I, O I, Na I, Mg II, 
Ca II, and Fe II. This may occur even if helium is present with a total 
mass of ~ 1 Mq , and thus comparable to the models characterised 
by a surface helium mass fraction of ^85 %. Early-time spectra are 
thus an important tool to determine whether the progenitor star is a 
WN (in a binary system) or a WCAVO (single) star. 

(viii) Our simulations of W-R-star explosions at low metallicity 
show weaker metal lines, and in particular those from Fe II. Be- 
cause spectra of SNe Ilb/Ib/Ic taken prior to or at the time of re- 
brightening are unaffected by the explosively nucleo-synthesised 
metals, quantitative spectroscopy of early-time SN observations 
may offer a powerful and accurate means to directly infer the 
metallicity of the primordial gas from which the progenitor of 
the SNt formed. This would also be a direct way of constrain- 
ing the effect of metallicity on stellar evolution, by passing alter- 
nate measurements performed on the galaxy env i ronment of the 
SN llKewlev&Ellisorj[2003 : iModiaz et all (20081 : iLevesque et all 
I2OIOI) . 

(ix) SN classification suffers from a great bias due to the time 
of observation. Many events classified as lb could have been lib's 
if discovered earlier, when Ha appears as a strong P-Cygni profile. 
Most SNe Ilb/Ib/Ic are detected during the brightening to peak, or 
sometimes so late that even the peak is missed. Failing early detec- 
tion, seeking He 1 10830 A is key to address the presence of Helium 
in the progenitor star, and discriminate between a type lb and a type 
Ic SN. We are cuiTently in an embarrassing situation, with SNe lb 
(Ic) showing signs of hydrogen (helium). 

(x) Our set of binary models, systematically characterised by 
low-mass ejecta, ^®Ni present out to ~2500km s^^, and a very 
large helium surface abundance, seem to reproduce the key prop- 
erties of SNe lib and lb. Our single star models, which give rise 
to SNe Ic owing to the lack of helium lines in their spectra, 
have lower but still significant (X (He)<50%) surface He mass- 
fractions. Not included in our sample here, the lowest-mass mas- 
sive stars in binary systems are expected to produce pre-SN ob- 
jects with very little surfac e helium, likely producing Type Ic SNe 
dWellstein & Langei|[l999l) . In this context, SNe lib and lb could 
stem from binary-star evolution, and SNe Ic could stem from both 
single- and binary-star evolution with an obvious bimodal d istribu- 
tion in ejecta mass (the recent analysis of lOrout et al .I2OIOI reflects 
this in part). In a forthcoming study, we will include such progeni- 
tors, as well as allow for the presence of ^''Ni in all ejecta to com- 
pare the light-curve morphology for all resulting SNe. 

(xi) Owing to mass loss through mass transfer to a companion 
star, massive stars evolving in a binary system can become W- 
R stars at a much lower main-sequence mass, p erhaps as low as 
IOMq jVanbeveren et al.lll998l : lYoon et al.ll201(3l) . A single star of 
the same mass would die as a RSG (it would also give rise to a 
Type II-P SN). Because the initial-mass function favors the forma- 
tion of lower-mass stars, we expect a large number of low-mass 
low-luminosity W-R stars. Paradoxically, most W-R stars are in- 
feiTed to stem from stars with a main-sequence mass in excess 
of 25 Mq , p otentially extending to huge masses of a few 100 Mq 
jCrowther et al. 2010), generically characterised by a large lumi- 
nosity of lO^-lO '' Lq and large mass-loss rates of ~10^''Mq yr^^ 
( ICrowthetll2007h . 

One striking example of this mismatch is the hydrogen-rich WN 
progenitor star corresponding to model Bmil8mf3p79zl. It has a 
luminosity of ~7OOOOL0, a surface temperature of 29000 K, an 



© 2010 RAS, MNRAS OOQ.fTltlB] 



Radiative transfer ofSN Ilb/Ib/Ic ejecta 2 1 



escape velocity of 380km s^^, is at 57% of the Eddington lumi- 
nosity, and is expected to have a feeble wind mass loss rate of 
~5xlO~ M0yr~^. In contrast, known WNh stars are amongst 
the most luminous and most massive W-R stars observed today. 
Hence, the higher-mass higher-luminosity W-R stars we know to- 
day seem to have little in common with the progenitors of cuiTently 
observed SNe Ib/c. This issue is in fact not incompatible with ob- 
servations. Low-luminosity low-mass W-R stars in binary systems 
are sitting next to a bright and more massive companion, to which 
the primary's envelope was largely transferred. The W-R luminos- 
ity is dwarfed by that of the OB-star companion and is consequently 
not easily identifiable. Their longer-term descendan ts, the Be/X-ray 
binar i es, are substantial e vidence of this scenario dWellstein et al.l 
l200ll : |Petrovic et al.ll2005l) . 

(xii) If the majority of SNe Ilb/Ib/Ic progenitors are from low- 
mass massive stars in binary systems, the notion that the sequence 
II-P — i- II-L — >■ lib — >■ lb — >■ Ic represents one of increasing main- 
sequence mass and increasing mass loss needs revision. Instead, 
one could understand this trend merely by invoking single- or 
binary-star evolution. Namely, stars with a main-sequence mass 
<2OM0 would make SNe II-P and II-L if single, and SNe Ilb/Ib/Ic 
if part of a binary system (this does not exclude that a subset of SNe 
Ic could come from single massive stars). Numerous studies have 
focused on understanding how such massive objects could lose suf- 
ficient amount of mass and produce a successful explosion despite 
their highly-bound high-density cores. If most SNe Ilb/Ib/Ic stem 
from lower-mass massive stars, these issues are largely irrelevant. 
In this context, high-mass W-R stars would collapse to a black hole 
without any associated SN or GRB signal. 

(xiii) The possibility that most, if not all, SNe Ilb/Ib/Ic progen- 
itors stem from low-mass massive stars in binary systems would 
modify the scene for the core-collapse explosion mechanism. Such 
low-mass progenitors have small and loosely-bound cores that will 
be no more difficult to explode than those of SNe II-P progenitors. 
In contrast, if SNe Ilb/Ib/Ic came from high-mass massive stars, 
a very powerful explosion mechanism, which is presently lacking, 
would be needed to successfully eject their massive highly-bound 
envelopes. If all core-collapse SNe come from low-mass massive 
stars, either single of binary, the challenge for the explosion mecha- 
nism is to be viable for 10-20 M© main-sequence mass stars, some- 
thing that may be more easily attained. Such SNe would also expe- 
rience little or no fallback, corroborating the observation that low- 
luminosity ''''Ni-p oor SNe II-P come from 8-IOMa massive stars 
(e.g., SN 2005CS. iMaund et alj|2003 ; iLi et alj[20o3 ; SN 2008bk, 
iMattila et al1l2008ir Such binary-star progenitors would thus even- 
tually make an important contribution toward the oxygen enrich- 
ment of the interstellar medium. 

Despite the recent observational efforts, our view of SNe 
Ilb/Ib/Ic (and la) is compromised by the lack of early-time spec- 
troscopic and photometric observations, many days prior to peak 
when the SN is faint and the photosphere still resides in the outer 
layers of the progenitor star. Such observations would offer useful 
constraints on the progenitor surface composition, the metallicity 
of the progenitor molecular cloud, and consolidate the determina- 
tion of the ejecta mass. Forthcoming blind deep full-sky survey will 
allow us to resolve these biases, by capturing a near-complete sam- 
ple of stellar explosions covering from a very early post-explosion 
time up to late times. 

Non-LTE time-dependent radiative-transfer modelling of the 
type presented here allows the simultaneous computation of spectra 
and light curves with the same high level of physical consistency. In 



the future, with the added treatment of heating/excitation/ionisation 
from non-thermal electrons, this approach will permit a better de- 
termination of the level of ''^Ni-mixing in SN ejecta and a more 
robust modelling of SNe Ilb/Ib/Ic, concerning both the progenitor 
as well as the explosion properties. 
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